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treatment  of'  these  topics,  has  been  presented  in  an  earlier  report.  Therein  it  was 
shown  that  the  surface  characteristics  of  the  scatteror  can  be  reluted  to  the  en~ 
semble  average  power  spectral  density  of  the  scattered  radiation,  Huwcver,  since 
the  measurement  of  the  power  spectral  density  is  difficult  and  a  more  practical  ob¬ 
servable  is  the  ensemble  average  scattered  intensity  distribution,  an  analytical 
expression  for  the  latter  is  developed  from  basic  principles  and  related  to  some 
average  properties  of  the  scattering  surface. 

By  assuming  a  Gaussian  distribution  ior  the  surface  elevation,  the  average  intensity 
distribution  is  related  to  the  second-order  statistics  of  the  surface.  It  is  shown 
that  this  relationship  is  similar  to  that  obtained  in  the  case  of  high-frequency 
microwave  backscatter.  The  mathematical  parameters  related  to  the  covariance  func¬ 
tion  of  the  surface  are  shown  to  lead  to  the  mean  square  elevation,  slope,  and 
curvature  of  the  surface. 

In  order  to  apply  these  theoretical  results  to  the  ocean  surface,  n  wind-dependent 
two-dimensional  power  spectral  density  for  the  ocean  surface  elevation  (based  on  a 
large  amount  of  empirical  data)  is  assumed.  The  surface  covariance  function  is  de¬ 
termined  and  it  is  shown  that  the  average  scattered  intensity  distribution  is  de¬ 
termined  solely  by  the  slope  statistics  of  the  surface  for  the  case  of  optical 
illumination.  An  analytical  expression  for  the  scattered  intensity  distribution  is 
derived  for  a  general  geometry.  The  detailed  distribution  is  calculated  for  the 
backscatter  geometry  and  its  behavior  is  examined  as  a  function  of  wind  speed.  It 
is  found  that  as  the  wind  increases,  the  bochscattcred  Intensity  is  reduced  and 
spread  over  lurger  angles.  It  is  inferred  that  the  downwind  slopes  increase  only 
slightly  faster  than  the  crosswind  slopes, 

A  reasonable  parametri<M-d  form  is  assumed  for  the  spectral  shape  of  a  perturbation. 

A  parametric  study  was  porforisc-d  and  it  was  found  that  the  effect  of  the  perturba¬ 
tion  on  the  intensi  distribution  is  most  pronounced  in  the  capillary  region.  The 
energy  content  of  tho  perturbation  required  to  see  significant  changes  is  on  the 
same  order  as  tho  energy  content  of  the  natural  spectrum  in  the  capillary  region 
and  insignificant  compared  to  the  total  energy  of  the  waves.  Tho  viscous  cutoff 
does  not  havo  much  effect  on  the  intensity.  The  angular  shipe  and  the  orientation 
of  the  perturbation  affect  mainly  tho  symmetry  of  the  5 n tens ttv  pattern  with  respect 
to  the  wind  direction.  Some  mathematical  details  regarding  tho  surface  covariance 
function  arc  presented  in  tho  two  appendices. 
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SUMMARY 


There  has  been  considerable  controversy  in  the  optical  literature 
concerning  the  theory  of  ’’speckle  patterns"  resulting  from  the  diffuse 
scattering  of  coherent  light  and  the  information  about  the  scatterer 
obtainable  from  the  speckle  pattern.  The  resolution  and  reconciliation 
of  these  theoretical  differences  and  the  application  of  the  theory  to 
scattering  from  the  ocean  surface  were  the  objectives  of  this  study. 

A  brief  historical  review  of  some  previous  theoretical  work  and  a 
description  of  the  apparent  discrepancies  which  have  now  been  reconciled 
are  given  (Section  II).  A  more  complete  treatment  of  these  topics,  has 
been  presented  in  an  earlier  report.  Therein  it  was  shown  that  the  sur¬ 
face  characteristics  of  the  scatterer  can  be  related  to  the  ensemble 
average  power  spectral  density  of  the  scattered  radiation.  However, 
since  the  measurement  of  the  power  spectral  density  is  difficult  and  a 
more  practical  observable  is  the  ensemble  average  scattered  intensity 
distribution,  an  analytical  expression  for  the  latter  is  developed  from 
basic  principles  and  related  to  some  average  properties  of  the  scattering 
surface  (Section  II-C). 

By  assuming  a  Gaussian  distribution  for  the  surface  elevation,  the 
average  intensity  distribution  is  related  to  the  second-order  statistics 
of  the  surface.  It  is  shown  that  this  relationship  is  similar  to  that 
obtained  in  the  r>rc  of  high-frequency  microwave  backscatter  (Section 
II-D).  The  mathematical  parameters  related  to  the  covariance  function 
of  the  surface  are  shown  to  lead  to  the  mean  square  elevation,  slope, 
and  curvature  of  the  surface  (Section  II-E), 
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In  order  to  apply  these  theoretical  results  to  the  ocean  surface, 
a  wind-dependent  two-dimensional  power  spectral  density  for  the  ocean 
surface  elevation  (based  on  a  large  amount  of  empirical  data)  is  assumed 
(Section  III-B).  The  surface  covariance  function  is  determined  and  it  is 
shown  that  the  average  scattered  intensity  distribution  is  determined 
solely  by  the  slope  statistics  of  the  surface  for  the  case  of  optical 
illumination  (Section  III-C).  An  analytical  expression  for  the  scattered 
intensity  distribution  is  derived  for  a  general  geometry.  The  detailed 
distribution  is  calculated  for  the  backseat tor  geometry  and  Its  behavior 
is  examined  as  a  function  of  wind  speed  (Section  IIX-D).  It  is  found 
that  as  the  wind  increases,  the  backscattered  Intensity  is  reduced  and 
spread  over  larger  angles.  It  is  Inferred  that  the  downwind  slopes  in¬ 
crease  only  slightly  faster  than  the  crosswind  slopes. 

A  reasonable  parametrized  form  is  assumed  for  the  spectral  shape  of 
a  perturbation  (Section  III-E).  A  parametric  study  was  performed  and  it 
was  found  that  the  effect  of  the  perturbation  on  the  intensity  distribution 
is  most  pronounced  in  the  capillary  region.  The  energy  content  of  the 
perturbation  required  to  see  significant  changes  is  on  the  same  order  as 
the  energy  content  of  the  natural  spectrum  in  the  capillary  region  and 
insignificant  compared  to  the  total  energy  of  the  waves.  The  viscous 
cutoff  does  not  have  much  effect  on  the  intensity.  The  angular  shape 
and  the  orientation  of  the  perturbation  affect  mainly  the  symmetry  of 
the  intensity  pattern  with  respect  to  the  wind  direction  (Section  II1-F), 
Some  mathematical  details  regarding  the  surface  covariance  function  are 
presented  in  the  two  appendices. 
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ISE  2618,  The  technical  direction  of  the  program  was  provided  by 
Dr.  Richard  F.  Hogiund,  Chief,  Advanced  Concepts  Division,  Strategic 
Technology  Office,  Defense  Advanced  Research  Projects  Agency,  Arlington, 
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I  INTRODUCTION 


A.  Motivation 

Scalar  diffraction  theory  has  been  applied  with  considerable  quanti¬ 
tative  success  to  numerous  microwave  and  optical  problems  for  which 
spatial  coherence  is  preserved.  However,  confidence  in  theoretical  pre¬ 
dictions  is  not  high  for  problems  involving  scattering  from  rough  sur¬ 
faces.  Often,  analytical  solutions  to  scattering  problems  can  be  obtained 
only  by  invoking  simplifying  approximations  that  are  known  to  be  rather 
crude.  Moreover,  experimental  testing  of  scattering  theory  is  difficult 
because  of  the  dual  problem  of  not  only  obtaining  a  surface  for  which 
the  theoretical  approximations  are  valid  but  also  specifying  the  statis¬ 
tical  characteristics  of  the  surface.  In  the  case  of  optical  scattering 
from  a  rather  calm  sea,  the  usual  gentle  slope  and  gentle  curvature 
assumptions  seem  to  be  reasonably  well  justified  and  one  would  expect 
quantitative  predictions  to  bo  moderately  accurate  provided  the  Fresnel 
and  Fraunhofer  conditions— as  well  as  the  small  angle  condi tions — are 
met.  Although  experimental  testing  of  these  predictions  cannot  yet  be 
accomplished  because  there  is  no  technique  now  available  to  measure  the 
statistical  characteristics  of  the  sea  with  the  accuracy  and  resolution 
that  is  required,  such  techniques  exist  in  principle,  and  when  available 
in  hardware  form  could  provide  the  means  to  perform  this  important  test. 
The  recognition  of  these  facts,  and  the  desire  to  develop  practical 
methods  of  utilizing  the  powerful  coherent  light  sources  that  are  now 
available  for  measuring  ocean  wave  statistics,  have  motivated  this 
reseurch. 
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\  B.  Outline  of  the  Report 
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There  has  been  considerable  controversy  in  the  optical  literature 
concerning  the  theory  of  "speckle  patterns"  from  diffusely  scattered 
coherent  light.  A  brief  historical  review  of  some  previous  theoretical 
work  is  included  here  as  well  as  a  description  of  some  of  t lie  apparent 
discrepancies  which  now  have  been  reconciled.  This  is  done  in  order  to 
demonstrate  the  foundation  of  the  theoretical  approach  taken  here. 

Starting  with  the  Huy gens- Fresnel  principle,  a  principle  familiar 
to  optical  physicists,  we  develop  an  analytical  expression  for  the  ob¬ 
served  ensemble  average  scattered  Intensity,  and  relate  it  to  statistical 
characteristics  o»  the  scattering  surface.  This  result  is  compared  with 
the  microwave  backseat  ter  cross  section  that  is  derived  from  the  Helmholtz 
integral  formalism  and  leads  to  a  new  physical  Interpretation  of  the  bnek- 
sentter  cross  section.  We  place  a  physical  interpretation  on  the  im¬ 
portant  mathematical  parameters  that  determine  the  scattered  intensity. 

Finally,  because  sufficiently  accurate  experlmcn.nl  data  describing 
the  statistical  characteristics  of  the  ocean  are  not  available,  we  use  a 
recently  reported  theoretical  description  of  the  ocean-wave-height 
power  spectral  density  to  make  numerical  calculations  of  the  scattered 
intensity  distribution  as  a  function  o)  wind  speed  and  direction.  In 
addition,  we  calculate  the  scattered  intensity  distribution  when  certain 
assumed  perturbations  having  arbitrary  angular  orientation  with  respect 
to  the  wind  direction  are  superimposed. 


References  are  listed  at  the  end  of  this  report, 


C.  Definition  of  the  Problem 

The  problem  that  is  addressed  in  this  report  can  be  outlined  with 
the  help  of  Figure  1.  A  plane  monochromatic  wave,  K  (u,v,z),  is  incident 
on  the  surface  of  the  ocean  at  the  origin  of  the  coordinate  system, 
(u,v,z).  The  mean  plane  of  the  ocean  surface  is  taken  to  be  the  (u,v) 
plane,  and  the  angles  A  and  B  define  the  plane  of  incidence  and  angle  of 
incidence  respectively.  The  observation  plane,  the  (x,y)  plane,  is 
parallel  to  the  <u,v)  plane  and  is  located  a  distance  h  above  it.  The 
observation  point,  Q  ,  is  defined  either  by  the  coordinates  (x,y,h),  or 
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the  vector  r^,  and  angles  a  and  0.  The  incident  plane  wave  has  the 
specific  form 

=  g(u,v,z)EQ  exp^j  k^  •  r^  (1.1) 

where  g(u,v,z)  is  an  amplitude  aperture  function  limiting  the  extent  of 
the  plane  wave  and 


“ktsin  B  cos  A  i  +  sin  B  sin  A  3  +  cos  B  k] 


-ui  +  v3  +  zk 


(1.2) 

(1.3) 


with  (i,j,k)  being  unit  vectors  along  the  (u,v,z)  directions  respectively. 


The  arbitrary  aperture  function  g(u,v,o)  defines  the  finite  portion 
of  the  ocean  that  is  illuminated  and  permits  infinite  limits  to  be  used 
on  integrals  over  the  (u,v)  plane.  For  example,  the  power,  P^,  incident 
on  the  (u,v)  plane  is  given  by 


(1.4) 


where  Z  is  the  Impedance  of  free  space  (120  n  ohms).  In  the  specific 
o 

case  for  which  the  incident  beam  has  a  Gaussian  cross  section,  the  Inte- 

2  2 

gral  Eq.  (1,4)  has  the  value  rrq where  Eq  is  the  peak  intensity  and  q 
is  the  e  ^  intensity  radius. 


Figure  2  will  aid  in  further  defining  the  geometry  of  the  scattering 
problem.  The  point  is  located  on  the  surface  of  the  ocean  within  the 

aperture  function  g(u,v,o)  and  has  coordinates  [u,v.w(u,vt , .  where  w(u.vl 
is  the  elevation  of  the  water  surface  above  the  mean  (u,v)  plane.  Alter¬ 
natively,  the  point-  Q,  may  be  defined  by  the  vector  r  and  the  angles  a 
and  b. 
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II  GENERAL  THEORY 

A.  Historical  Background  and  Review 

After  the  operation  of  the  first  CW  He-Ne  laser,  considerable 
interest  arose  in  the  speckled  appearance  of  typical  diffuse  surfaces 
illuminated  by  coherent  light.  It  was  soon  recognized  that  the  origin 
of  these  granular  patterns  lay  in  the  optical  roughness  of  the  surfaces 
from  which  the  light  was  scattered.2* 3 

The  first  detailed  statistical  treatment  of  this  phenomenon,  now 
known  as  laser  speckle,  was  that  of  Goodman/'  Treating  the  scattering 
or  reflecting  surface  as  a  statistical  ensemble  of  many  point-scatterers 
with  independent  and  uniformly  distributed  phases  on  (0,  2 n),  Goodman 
showed,  among  other  results,  that  the  power  spectral  density  (or  Wiener 
spectrum)  of  the  intensity  pattern  of  reflected  light  was,  up  to  scaling 
factors,  identical  with  the  autocorrelation  function  of  brightness  (i.e., 
intensity)  distribution  across  tho  scattering  spot. 

Goldfischer  independently  derived  the  same  result  using  essentially 
the  same  model  as  Goodman.4  Again  the  scattcrers  were  assumed  to  be 
infinitesimal  in  size  and  independent,  and  were  taken  to  have  phases 
uniformly  distributed  on  (0,2n).  This  assumption  is  mathematically 
equivalent  to  assuming  that  the  spatial  coherence  of  the  fields  behaves 
similarly  to  a  ^-function;  that  is,  the  spatial  coherence  has  value  unity 
for  zero  displacement  and  value  zero  for  finite  displace*  ant. 

I rt  an  extension  of  his  previous  results,  Goodman3  later  discussed 
the  properties  of  speckle  patterns  as  they  affect  the  performance  of 
optical  radars.  In  an  appendix  to  this  paper,  he  pi'esented  a  generalized 


7 


theory  that  allowed  n  finite  correlation  area  to  exist  in  the  fields  at 
the  scattering  surface.  Thus,  the  assumption  oi  infinitesimal  scattcrers 
was  removed,  but,  significantly,  the  phase  of  the  scattered  light  at  any 
point  on  the  object  was  still  assumed  to  be  uniformly  distributed  on 
(0,2o).  Hence,  the  theory  applies  only  for  surfaces  that  are  quite  rough 
on  the  scale  of  the  wavelength  of  the  illumination. 

Following  this  early  work  that  dealt  with  the  Intensity  patterns 
existing  in  the  scattered  light,  other  investigators — notably  Lowenthal 
and  Arsenault" — have  studied  the  properties  of  similar  patterns  observed 
in  the  images  of  coherently  illuminated  diffuse  objects.  Those  treat¬ 
ments  have  assumed  that  the  complex  fields  in  the  image  plane  obey  cir¬ 
cular  Gaussian  statistics;  i.e, ,  the  real  and  imaginary  parts  of  the 
field  are  zero  mean,  uncorrelated  Gaussian  random  variables,  with  identical 
variances.  Such  an  assumption  is  valid  if  the  correlation  area  of  the 
diffusely  transmitted  or  collected  wavefront  is  small  compared  with  the 
area  of  the  point-spread  function  of  the  imaging  system  and  if  the  phase 
of  the  light  at  each  point  on  the  object  is  approximately  uniformly  dis¬ 
tributed  on  (0,2n).  Since  the  current  discussion  will  be  limited  to 
speckle  patterns  observed  by  direct  measurement  of  intensity  in  the  scat¬ 
tered  light,  without  any  intervening  optics,  the  problem  studied  by 
Lowenthnl  and  Arsenault  will  not  be  addressed. 

The  next  important  piece  of  work  dealing  with  directly  observed 
speckle  patterns  is  that  of  Crane, 5  Hu  derives  results  that  appear  in 
various  respects  to  conflict  with  earlier  results  ol  Goodman'1’ n  and 
Goldfischer.  Crane’s  comments  are  directed  at  the  paper  by  Guldfischer, 
but  his  chief  criticism  is  as  follows:  Goldfischer  assumes  that  the 
field  produced  by  a  small  subsection  of  a  diffuse  surface  is  proportional 

t r\  |  m n  »» a  rnrv  t  nnha  n  j*  l  ]m  t  utiljcon  i  4  nn  •  hit  |  fni*  r«nUnti|vn  I 

light,  the  field  must  be  directly  proportional  to  the  area  because  the 
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field  contributions  odd  on  an  amplitude  basis.  Thus,  Crane  concludes 
that  Goldfischer  used  a  physically  incorrect  model  and  that  this  is  the 
reason  for  the  apparent  differences  between  the  results  of  the  two 
theories.  It  should  bo  mentioned  that  Crane's  theory  has  the  merit  of 
allowing  for  a  finite  correlation  area  of  the  diffuse  surface.  Of  the 
work  preceding  Crane’s,  only  that  of  Goodman'1  allowed  for  this  possibility. 

Crane's  criticism  of  Goldfischer' s  work  was  answered  in  a  letter  by 
Arsenault.2  Arsenault  asserts  that  the  differences  between  Goldfischer’ s 
and  Crane’s  results  arise  because  the  former  theory  derives  statistical 
average  results  while  the  latter  is  entirely  a  deterministic  theory. 
Arsenault  addresses  his  comments  only  to  the  early  part  of  Crane's  paper 
and  in  particular,  to  that  portion  through  Crane's  Eq.  (IS).  In  a  later 
portion  of  his  paper  leading  to  his  Eqs.  (38)  and  (39),  Crane  has  derived 
expressions  for  the  statistical  average  power  spectral  density,  the  same 
quantity  derived  by  Goldfischer,  The  important  differences  between  these 
two  expressions  wore  not  reconciled  by  Arsenault's  comments.  However, 
in  1973  Goodman40  showed  that  the  discrepancies  are  only  apparent  and  not 
real,  and  that  they  have  resulted  from  improper  interpretation  of  results 
rather  than  improper  development  of  the  theory.  Specifically,  Crane's 
result  as  expressed  by  his  Fq.  (39)  can  only  be  applied  to  relatively 
smooth  surfaces,  that  is,  surfaces  for  which  the  rms  surface  roughness 
is  much  less  than  one  wavelength  of  the  scattered  radiation.  On  the 
oth^r  hand,  the  results  of  Goodman's  and  Goldfischer’ s  analysis  can  only 
be  applied  to  truly  rough  surfaces,  that  is,  to  surfaces  for  which  the 
rms  surface  roughness  is  much  greater  than  one  wavelength  of  the  scattered 
radiation.  Moreover,  Goodman40— starting  with  Crane's  Eq,  (38) — derived 
an  approximate  expression  for  the  statistical  average  power  spectral 
uensiiy,  valid  for  truly  rOu^ii  surfaces,  that  is  entirely  consistent 
with  all  previous  theoretical  work.  Thus,  the  apparent  discrepancies 
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have  been  reconciled  and  all  theories  are  mutually  consistent  when 
properly  applied  to  their  regions  of  validity. 


B.  Approximations  and  Assumptions 

The  mathematical  statement  of  the  Huygens- Fresnel  principle  is 


E(x.y) 


_  _1_ 
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Implicit  in  the  use  of  this  principle  are  approximations  and  assumptions 
similar  to  those  used  if.  the  Helmholtz  integral  formalism — the  usual 
starting  point  for  microwave  scattering  calculations.  A  more  complete 
discussion  of  these  approximations  and  assumptions  will  be  found  in 
Beckmann  and  Spizzichino.  *  - 

As  in  the  case  of  ail  scalar  theories,  the  results  obtained  here 

are  valid  only  for  small  angles  of  incidence  and  small  bistatic  angles. 

We  make  the  approximation  that  the  field  E  (u,v,w)  reflected  from  the 

surface  can  he  represented  by  the  product  u, v)  E  (u^VjW),  where  Tj(u,v) 

is  the  surface  reflectance.  This  is  a  reasonable  approximation  since 

the  minimum  water  wavelength  A  is  known  to  be  much  larger  than  the 

min 

optical  wavelength  The  amplitude  reflectance,  l)(u,v),  is  taken  to 
be  -|T|(B)  |,  where  ^(B)  is  the  Fresnel  amplitude  reflectance  at  the  inci¬ 
dence  angle  B;  the  surface  elemental  area,  dS  ,  is  taken  to  be  the  same 

S' 

as  the  elemental  area  dudv  in  the  (u,v)  plane.  The  obliquity  factor  is 
taken  to  be  cos  g.  The  preceding  three  approximations  are  reasonable 
for  a  rath  calm  sea  with  gentle  slopes.  Finally,  the  factor  Jr^  -  r^| 
is  approximated  by  r  in  the  denominator  of  Eq.  (2,1),  and  also  by  the 
first  few  terms  of  its  binomial  expansion 


2  2 

i  i  xti  1  yv  u  I-  v 
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(2.2) 


in  the  exponent  of  Eq.  (2.1).  In  this  regard  it  should  be  noted  that 
there  are  different  ways  of  grouping  the  terms  in  the  binomial  expansion 
that  result  in  slightly  different  approximations.  In  one  case,  the  small 
angle  is  approximated  by  its  tangent  while  in  the  other  case  it  is  approxi¬ 
mated  by  its  sine.  The  Huygens- Fresnel  principle,  in  its  approximate  form, 
is  stated  here  as  Eq.  (2.3). 
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The  basic  statistical  assumption  is  made  here  that  the  ocean  wave 
height  statistics  are  wide-sense  stationary,  at  least  over  the  patch  of 
ocean  illuminated,  and  the  wave  heights  are  Gaussian  distributed, 1 3 
The  effects  associated  with  sampling  only  a  finite  patch  of  the  ocean 
both  in  the  wave  height  statistics  and  the  optical  field  statistics  are 
accounted  for.  The  statistics  of  the  scattered  fields  at  the  observation 
point  are  taken  to  be  circular  complex  Gaussian.  Such  statistics  arise 
only  when  the  scattered  field  components  generate  a  complex  random  walk, 
with  approximate  uniform  phase  associated  with  the  various  field  contribu¬ 
tions.  This  amounts  to  assuming  that  the  rms  surface  roughness  is  com¬ 
parable  with  or  greater  than  the  wavelength  of  the  scattered  radiation, 
an  assumption  that  is  well  justified  oven  for  a  calm  ocean.  It  does  mean, 
however,  that  the  present  results  cannot  be  applied  to  perfectly  smooth 
surfaces,  such  as  polished  mirrors. 
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C .  Ensemble  Average  Scattered  Intensity 

In  this  section  an  express  ion  is  derived  for  the  ensemble  average 
scattered  intensity  distribution,  (I(x,y)),  as  a  function  of  the  nor¬ 
malized  ensemble  average  spatial  coherence  function  of  the  fields  in  the 
(u,v)  plane.  Subsequently  it  is  shown  how  the  spatial  coherence  can  be 
expressed  in  terms  of  physically  observable  ensemble  average  surface 
characteristics. 

An  explicit  expression  for  the  field  E^(u,v,w)  can  be  written  with 
the  help  of  Eqs.  (1.1),  (1.2),  and  (1.3). 


E  (u.v.w) 
s  ' 


E(u,v)  exp 


-jk[*o(u,v 


)  t  w  cos  B 


(2.5) 


with 


i)o(u,v)  =  u  sin  B  cos  A  *■  v  sin  B  sin  A  *  ~  ^ux^  F  vy^j  (2. 


6) 


and 


E(u,v)  ^  -|1|(B)  jg(u,v,o)  Eq 


(2.7) 


Equation  (2.3)  now  may  be  restated  as 


E(x,y)  = 


2  Jkro  " 
cos  8  e 


|  E(u,v) 
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exp  |jk|jlf(u,v)  -  >io(u,v)  -  w(cos  B  t  cos  g)j  |  dudv 


(2.8) 


As  on  aid  to  interpreting  these  results,  note  that  the  Huygens- 
Fresnel  principle  may  be  stated  in  the  equivalent  form6 
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_  j  fcr  “ 

2  J  off 

E(x,y)  =  — - -  II  E  (u,v)  exp[jk$  (u,v>]  dudv  (2.9) 

jXh  I*  P  ° 


where  E^(u,v)  is  the  mathematically  equivalent  field  in  the  (u,v)  plane 
as  opposed  to  the  field  on  the  scattering  surface.  The  terminology 
"mathematically  equivalent"  field  is  used  here  to  emphasize  t.he  fact  that 
Ej,(u,v)  is  not  the  Physical  field  in  the  (u,v)  plane  but,  rather,  it  is 
equivalent  to  a  physical  field  which,  if  present  in  the  (u,v)  plane, 
would  produce  the  field  E(x,y).  A  comparison  of  Eq.  (2.9)  •  th  Eq.  (2.8) 
yields  the  result 

E^(u,v)  =  E(u, v)  exp|-jk^o(u,v)  -  w(cos  B  +  cos  g)j  j  .  (2.10) 

The  ensemble  average  scattered  intensity,  (I(x,y)),  is,  by  definition 
<I(x,y)>  k  <E(x,y)E*(x,y>>  (2,11) 


where  the  carets  are  used  to  indicate  that  an  ensemble  average  is  to  be 
performed.  Equation  (2.9)  may  be  used  in  Eq.  (2.11)  to  yield 

a<X,y>>  =7^lt[f  <Ep(VVi)Ep(vV2)> 

(\h)  JJJJ 

~C3 

x  sxpjjkfs^.vJ  -  j(a2,vJ]{  1u1<l»1du2<.v2  .  (2.12) 


The  quantity  in  carets  in  the  integrand  will  be  recognized  as  the  mutual 
intensity  J(u  ,v  ,u  ,v  )  of  the  equivalent  fields  in  the  (n,v)  plane,  and 

X  X  £  £ 

from  Eq.  (2.10)  is  given  by 

£(vvvv2)>  -  e(vviMvv2)  ■  -o(vva)]i 

^xpjjk^w^  -  w2j(cos  B  +•  cos  g)j^>  .  (2,13) 
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knowledge  and  the  approximation 
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Eq.  (2.18)  now  can  be  written  as 
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c,d 

and  the  notation  3  is  used  to  indicate  a  two-dimensional  Fourier  trans- 
a,b 

form  with  respect  to  spatial  frequencies  (a,b)  and  spatial  coordinates 
(c,d). 


The  approximation  indicated  by  Eq.  (2.13)  is  reasonable  as  long  as 
the  angle  subtended  by  the  aperture  function  at  the  observation  point  is 
^.'.aller  than  about  one  degree. 

The  Integral  in  square  brackets  in  Eq.  (2.20)  may  be  evaluated  with 
the  help  of  Eqs.  (2.2L) ,  (2.18),  (2.7),  and  (1.4). 
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where  P  is  used  to  indicate  the  reflected  power.  Finally,  Eq.  (2.20) 
r 

may  be  rewritten  ns 
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There  is  an  interesting  interpretation  of  Eq.  (2.23)  that  is  con¬ 
sistent  with  the  understanding  of  the  propagation  of  spatial  coherence. 

For  the  case  of  an  infinite  aperture,  Eq.  (2.15)  reduces  to 

6/u  ,v  ;Au,Av)\  =  jl|(B)  |2E“(v(Au,Av)> 

In  this  form,  (y(Au,Av))  will  be  recognized  as  the  space  invariant, 
normalized,  ensemble  average  spatial  coherence  function.  Thus,  Eq. 

(2.23)  indicates  that  it  is  the  spatial  coherence  that  primarily  deter¬ 
mines  the  distribution  of  scattered  intensity.  It  is  therefore  important 
to  be  able  to  relate  the  spatial  coherence  to  surface  characteristics  in 
a  manner  that  will  enable  numerical  calculations  to  be  made. 

It  is  possible  to  derive  expressions  for  other  ensemble  average  ob¬ 
servables,  such  as  spatial  coherence  and  power  spectral  density.14  The 
ensemble  average  spatial  coherence  in  the  observation  plane  can  be  shown 
to  be  the  Fourier  transform  of  the  aperture  function  on  the  scattering 
surface  and  is  therefore  not  related  to  surface  characteristics.  More¬ 
over,  while  the  ensemble  average  power  spectral  density  of  the  scattered 
radiation  can  be  related  to  surface  characteristics,  the  experimental 
technique  of  measuring  this  quantity  would  be  intensity  interferometry 
and  therefore,  in  a  practical  sense,  it  would  be  more  difficult  to  measure 
than  the  ensemble  average  scattered  intensity  distribution.  Finally,  it 
should  be  mentioned  that  since  the  intensity  averaged  over  many  "speckles" 
is  equivalent  to  the  intensity  for  temporally  incoherent  light,  the  ex¬ 
pressions  derived  here  are  valid  for  incoherent  quasi-monochromatic  light 
sources,  such  as  filtered  sunlight. 
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D.  Spatial  Coherence  as  a  Function  of  Surface  Characteristics 


The  probability  distribution  function  p[w(u,v)]  for  the  surface 
heights  w(u,v)  must  be  known  in  order  to  perform  ensemble  averages.  An 
assumption  that  is  in  reasonable  agreement  with  experimental  evidence12 
is  that  the  surface  height  is  a  stationary  zero-moan  Gaussian  random 
variable  with  distribution 


plw(u,v)l 


,  /rr-x-1  {  [W(U,V>]2) 

=  (o  /2n)  exp j-  - ^ - 1 


(2.24) 


where 


2 

o  “ 


r 

2  §  o . 

^  [w(u,v)l  p[w(u,v)]  dw(u, v)  i  (Jw(u,  v)  ] “ J  ,  (2.25) 
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The  covariance  function  R(Au,Av)  and  the  normalized  covariance  func¬ 
tion  p(Au,Av)  for  homogeneous  (stationary)  processes  are  defined  by 


U(Au,Av)  i  ( w(u,v)w(u  +  Au,v  f  Av)) 


(2.26) 


and 
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The  joint  probability  distribution  function  p[w^,w  ]  is  given  by 
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where,  for  convenience,  he  following  notation  has  been  adopted: 


w  ^  w(u,v)  w  i  w(u  »  A u,v  i  Av)  . 
1  ~ 


(2.29) 
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The  quantity  in  ensemble  brackets  in  Eq.  (2.17)  can  be  interpreted 
as  the  characteristic  function  for  the  distribution  p{w  and  thus 

Eq.  (2.17)  can  be  written 

(vt&u.*Av))  =  exp[-Jk(Au  sin  11  cos  A  +  Av  sin  B  sin  A)] 

X  expj~k2(cos  B  i  cos  3)2[R(0,0)  -  R(Au,Av)]j 

*  0X"[-J  5  (iuxo  "  4V0] 

(2  2  ) 

X  expj-k  (cos  B  +  cos  0)  [R(0,0)  -  R(Au,Av)]j  .  (2.30) 

Through  Eq.  (2.30)  the  spatial  coherence  function  is  related  to  the 
surface  height  covariance  function  which,  with  the  known  (or  assumed) 
probability  distribution,  completely  characterizes  the  surface.  Finally, 
using  Eq.  (2.30)  in  Eq.  (2.23),  we  get 

’  fi!" f  (2Z°Pr) 

Au,Av 

X  ^exp|-k2(cos  B+cos  j3)2[R(0,0)-R(Au,Av)]J^ 

\\r  ^Jf\\r  R  J  (2.31) 

It  is  of  some  interest  to  compare  this  result  with  the  high  fre- 

o 

quency  microwave  backscatter  cross  section  c  .  In  our  notation  the  bank- 
scattered  (monostatlc)  intensity,  I,  is  given  by  Frost  and  Jackson:1^ 


2Z  Pj 

.  o  i  o 

1  =  ~  2  ° 
4nr 


(2.32) 
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5  r. 


t  » 


P  i 


h  *• 


l 


M 


i-  .* 


r  t 


e  * 


l  : 


i  v 


i  < 
*  t 


e  I 
I  3 


where 


o  2/  4n 

^  =  |TK3)  1 


2 

A  /  J 


cos(2kAu  sin  3) 


(22  ) 

X  expj-4k  cos  g[R(0,0)  -  R(Au,Av)]j  dAudAv 


(2.33) 


When  the  backscatter  conditions,  B  =  3,  A  =  0,  x  =  sin  g,  y  =  0,  are 


used  in  Eq.  (2.23),  the  result  is 

2 


cos 


>Av 


.  0^2Z  P  \  AU,i 

sin  9>0  )>=-7^  3  t(y(Au,Av))J  .  (2.34) 


kf 


sin  B 


K  ’ 


-.0 


We  may  construct  a  backscatter  cross  section  a'  using  Eq.  (2.32)  for 
the  definition  and  Eq.  (2.34)  for  the  backscattered  intensity 

Au,  Av 


0  s  co 


s2  B|TU9)!2|~j  [<y(Au,Av)>] 


(2.35) 


sin  B 


0 


\ 


Since  the  covariance  function  is  even  in  Au  (for  Av  =0),  the  Fourier 
cosine  transform  can  be  used  in  place  of  the  Fourier  transform,  to  ob~ 
tain  from  Eq.  (2.35) 


00 


2  2  /  4tt 

a'  =  cos  0  I  "HO)  |  |~r  cos(2kAu  sin  0) 

\>.2' 

—00 


(22  ) 

X  expj-4k  cos  0[R(O,O)  -  R(Au  Av)]j  dAudAv  (2.36) 


which  is  identical  with  Eq.  (2.33)  except  for  a  factor  of  cos  8.  This 
factor,  which  is  approximately  unity  for  small  angles,  arises  from  the 
obliquity  factor  in  the  Huygens-Fresnel  principle,  and  appears  to  be  the 
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only  difference  between  our  results  and  the  results  of  more  rigorous 
scalar  diffraction  theory. 

K,  The  Surface  Covariance  Function 

Let  w(u,v)  be  a  realization  of  the  random  process  representing  the 
elevation  of  the  ocean  surface.  Evidently  w(u,v)  is  a  real  zero-moan 
random  variable.  The  surface  covariance  function  R(Au,Av)  of  the  random 
surface  elevation  is  defined  by 

K(u, v,Au,Av)  4  (w(u,v)w(u  1  ^\\,v  t  ,\v))  .  (2.37) 

When  the  process  is  assumed  to  be  homogeneous  (stationary)  in  the  wide 
sense,  It  depends  only  on  the  displacements  (Au,&v)  ami  not  on  the  co¬ 
ordinates  ( u ,  v ) ;  i .  e . : 

(w(u,v)w(u  (  ,ju,v  *  Lv))  =  R(£u,£v)  .  (2.38) 

Further,  R  must  be  symmetric  about  zero  displacement;  i.e.: 

R(Au, Av)  =  R(-Au,-../)  .  (2.39) 

Or,  in  polar  coordinates 

R(Ar,e>  =  R(Ar,9  '  ")  .  <2.40) 

Since  the  measurement  of  the  surface  elevation  is  feasible  over  only  a 
finite  area  of  the  ocean,  it  implies  the  presence  of  an  aperture  g  (  «,v), 
sc  that  iho  «u^rnpf»  olovnt, ion  is  nivcn  by 

w  (u,v)  -  w(u,v)g  (u,v)  (2, 41) 
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where  g  (0,0)  =  1,  g  («,<»)  =  0,  and  [lim(n  -  m)g  (u,v)l  =  1,  so  that 
a  a  a 

w  (u,v)  =  w(u,v)  as  a  -  co,  if  W  (l  ,2  )  is  the  Fourier  transform  of 
g  g  u'  v 

w  (u,v), 
g 


CO 

fl 

• 

w  (u,v)  exp  - 

j(u  2  +  vi 

)  dudv 

* 

(2.42) 

■L 

g  L 

0 

\  U 

v/J 

The  ensemble  average  power  spectral  density  of  the  functions  w  (u,v)  is 

g 

defined  as 


;(v0 6 


(2,43) 


CO 

|j  |8a(u,v 


■J. 

)  dudv 


If  R  is  continuous,  then  by  the  use  of  the  Wiencr-Khintchine  theorem,  it 
follows  that 


CO 

R(Au,Av)  =  S*{K>K)  oxp[j(^u  h  Av\,)]  d*u 

(2n)  y 


d2  (2.44) 
v 


where 


00 

rr  r 

^u^v)  =  h<Au,Av)  expj^“J^Au£u 

U 


t  Avj?  H  dAudAv  (2.45) 


'U  ,1  )  =  lim  S*(t  ,i  ) 
\  u  v7  g\  u  v/ 


(2.46) 


and  S  (2^,2^)  is  known  as  the  power  spectral  density  of  the  random  process 
w(u,v).  It  may  also  be  shown  that 


S-Vfl  ,j£  ]\  =  R(Au,Av)  - — - - — -  A  R  .  (2.47) 

1  g\  u  v/J  PC  i  /  \i2  ^  .  g 

jj  |e  (u,v)l  dudv 
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If  the  covariance  of  the  aper lured 


functions  w  (u,v) 

g 


is  defined  as 


It  (ll,V,AU,Av)  4  <W  (  U  ,  V  )  W  (U 

g  ’  '  e  g 


AU,V 


(2,48) 


then 


It  (Au.Av) 
g 


JL» 


It  (u,v,Au,Av) 
it 


dudv 


(2,49) 


These  equations  express  the  result  that  the  covariance  function 
It(An,Av)  and  the  power  spectral  density  S*(Z  ,1^)  ol'  the  random  process 
w(u,v)  form  a  Fourier  transform  pair.  The  covariance  function  It  is  not 

7  b 

the  Fourier  transform  of  S *(jJ  .1  ).  However,  the  space  average  ol  It 

guv  8 

forms  a  Fourier  transform  pair  with  8*(2  ,  S.  ).  It  is  to  be  noted  that 

II  V 

It  is  not  homogeneous  and  that  (lim(a  -  «>)R  |  -  It(Au,Av).  Figure  3  shows 
g  & 

the  geometrical  relationship  between  the  several  domains  between  which 

the  transforms  are  defined. 


Frequently,  the  surface  covariance  function  cannot  be  obtained  fiom 
the  power  spectral  density  since  the  Fourier  transformation  indicated  in 
Kq.  (2.44)  is  difficult  to  perform.  For  this  reason,  as  well  as  to  ob¬ 
tain  a  more  direct  physical  insight  into  the  various  operations,  we  de¬ 
velop  K  in  a  Taylor  series  about  the  origin. 


R( Au,Av) 


(2.50) 


Because  of  the  symmetry  implied  by  Kq.  (2.39),  it  can  be  shown  that  in 
hq.  (2.50).  the  terms  lor  wlucli  n  is  odu  are  identically  zer u.  In  ulbci 
words,  all  odd  order  partial  derivatives  ol  R  are  zero  at  the  origin. 
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FIGURE  3  GEOMETRICAL  RELATIONSHIPS  BETWEEN  TRANSFORM  DOMAINS 
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^  w*^.V  •”*!(:  r+n- ^.-"t^iT*- ■  v  .-^t 5 5^rr,‘v v£7yfc  5*£?*3eSV 


Writing  Eq.  (2.44)  in  vector  notation,  and  using  the  definition 


S(£)  =  S*(£)/(2tt)‘ 


(2.57) 


we  get 


=  1. 


R(Ar)  =  J  S(0,)  exp(.j  Ar  •  A)  d£ 

plane 

Substituting  in  the  right  hand  side  of  L’q.  (2.53),  we  get 


R(Ar) 


m  ,  v  2n  f 
11-0  plane 


2n 

S(£)(Ar  *  £)  d£ 


(.0 


2m  12 


(2m  l-  2)1 


(2.58) 


S(£)(Ar 


. . 2m ‘2  / .  .  \  ,  . 

Ar  •  l)  exp  f  j  Pq  •  £ )  d_£  .  ( 


2.59) 


plane 

If  we  use  polar  coordinates  to  describe  both  R(Ar,0)  and  S(£,o),  we  have 


Ar  •  £  =  Ar  i  cos(0  -  *) 


and 


R(Ar,6) 


-E  ff 

»=°  *4  -In 


2n 

S(£,®)[£  cos ( 0  -  ©)]  £  d£do 


1[ir2”'2/(2»  .  2):]  f  f 

J0  -TT 


(-l)nH  1[Ar2m}2/(2m  *  2)j]  j  {  S(£,e)l£  cos(0  -  o)J2ml“ 


x  expfj  p  l  cos  9)  £  d£do 


(2.60) 
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Ff£5 i  ;«,^i*«w,<«^W’\"T>^^^"^«7sa'5TT5<’^7W!5a 


or 


m 


*4r'<»  V0^”  '  2 


,  rrH  1 


X)  7  Mo  0(P>Ar2mH2  (2.61a) 


n=0 


2m  -i-  21  2mi-2 


M_  (0) 
2n 


00 

■fj 


.  2n 


where 

■JO  v 

S(£,©)[£  cos(e  -  0)]““  A  d£d©  .  (2.61b) 

'O  '-T» 

It  is  thus  seen  that  the  various  terms  in  the  expansion  of  Eq.  (2.56) 
correspond  to  various  moments  of  the  power  spectral  density.  These 
moments  may  also  be  related  to  the  ensemble  average  properties  of  the 
random  process  w(u,v).  The  first  term  of  the  expansion,  or  the  zeroth 


moment,  is 

CO  TT 

S(£,©)  £  d£d©  ,  (2.62) 

'o  "-n 

using  Eqs.  (2.61a)  and  60).  From  Eq.  (2.38)  it  is  seen  that 


co  r 

•  i  ( 


R(0,G)  -  ^(w(u,v))2^ 


(2.83) 


Thus  the  zeroth  moment  gives  the  mean  square  elevation. 


The  second  moment  is  given  by 

CO  Tt 

I 

M 

2 

~0  '-TT 

and  the  second  term  of  the  expansion  is 


co  r 

■l 


3  2 

S(£,©)  £  cos  (0  -  *)  d£dc 


(2.64) 


1  ..  .  2  IT/  _  \2_,  .1 

2  m2  ar  *  i  i(&  *  'P;  kv£;j 


p=0 


26 


M2  may  also  be  expressed  as 


B!  ■  [(©  )  “  #  '  '((slfS))  sl"  #  ~  ‘  ((s)  )  ,ll!  9]  •  <*•»»> 


Thus  the  second  moment,  Mt>,  is  seen  to  be  equal  to  the  mean  square  of  the 
slope  component  in  the  direction  0. 

The  fourth  moment,  M  ,  may  be  similarly  related  to  the  partial  de- 
rivatives  of  v:(u,v),  which  are  related  to  the  curvature  of  the  surface. 
Tlius , 
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III  CALCULATION  OF  THE  SCATTERED  INTENSITY  DISTRIBUTION 

A.  General 

The  theoretical  relationships  developed  in  the  foregoing  are  applied 
in  this  chapter  to  the  scattering  of  optical  radiation  from  the  suriace 
of  the  ocean.  A  specific  power  spectral  density  is  assumed  for  the  sur¬ 
iace  elevation  of  the  natural  sea  surface,  and  an  expression  is  obtained 
for  the  scattered  intensity.  Subsequently,  a  specific  form  is  assumed 
for  the  power  spectral  density  of  a  perturbation  on  the  ocean  surface. 

An  expression  is  obtained  for  the  intensity  distribution  due  to  the  per¬ 
turbed  ocean  surface.  Numerical  calculations  are  made  to  evaluate  the 
effects  of  the  various  parameters  of  the  perturbation. 

B.  The  Power  Spectrum  of  Ocean  Surface  Elevation 

In  a  recent  report’  a  large  amount  of  experimental  data  were  digested 
to  develop  an  analytic  representation  of  the  spectrum  of  a  wind  roughened 
sea.  This  spectrum,  which  is  a  function  of  wind  velocity,  has  been  used 
in  the  present  study  as  the  power  spectral  density  of  the  surface  eleva¬ 
tion  of  ocean  waves. 

For  a  friction  velocity  u,  greater  than  a  minimum  value  u.  =  12  cm  s, 

*  *m 

the  spectrum  has  been  defined  by  the  following  equations. 


Masses ag 


S  *(£.«) 


(2r-) 


=  sa,c) 


S(£)F (2,*)  for 
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0  <  l  <  CO 


<?  O  <  TT 


(3.1) 


■l  (  r  2 ,  2  n  | 

-  (-1  3~)  cos  -r  1 1  -  expl-g  2 1  U  j 

*  ^  [l  ■  ^(u*)  cos  2o]  expf-^-A-1] 


(3.2) 


/  4\  r  2,  2  4 

^a/2 a  j  expj-Sg  i  U 

]  ° 

).  ■ 

?i 

(3.3a) 

/  1/2  7'2\ 

“'Pi  £  ) 

£! 

l  ■ 

*2 

(3. 31}) 

nl>  (m!>''-P) 

*2 

■  }, 

C3 

(3.3c) 

ai)/^2jf,*  j 

*3 

l 

a 

(3.3d) 

6  ,/  1(>\ 
aD£  /(2A  1 

l 

e  - 

0 

(3. 3e) 

,  =  0.74,  £  *  0.3588 

-1 

cm  , 

0.9437 

cm  1 ,  and 

U  -  the  wind  velocity  at  an  altitude  o'.  19.5  meters 

=  -2.5  u*|ln[(°.68‘,/u*)  ‘  4.28  x  10"5  u*  "  0.0443]  -  in  1950 1  (3.< 


*1 ■  "*y 


l ‘2  ,  1  '6 

i  =  0.5726  u  i.  I) 
m 


(3.5) 


(3.6) 


9.  =  («.,'-) 

m 

=  3.672  cm" 


0  -  ^1.27!  •  0, 


0268  ir  •  6,03  y  JO 
* 


(3.7) 


(3.8) 
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a,  ■=  0.214  log  u  0.036  .  (3.10) 

Thus  it  is  scon  that  S(£,<s)  is  dependent  on  the  1  fiction  velocity, 
or  indirectly  on  the  wind  speed. 

In  this  description  of  the  spectrum,  the  direction  ol  the  wind  has 
been  assumed  to  coincide  with  the  u-axis  or  the  direction  c  a  u  =  0. 

We  note  that 


j  F(£,c)  dc  ^  1  .  (3.11) 

In  Figure  *1,  l  S(£)  is  plotted  against  }..  The  area  under  the  curve 

is  equal  to  the  mean  square  elevation  according  to  Kq,  (2.62),  Figure  5 

presents  the  proven  spectral  density  on  a  log-log  plot,  accentuating  the 

dependence  on  wind  speed  in  regions  of  interest  to  this  study.  It  is  to 

be  noted  that  a  negligible  contribution  is  made  to  surface  roughness 

from  wave  numbers  larger  than  0.01  cm  Figure  6  shows  the  factor 
2  2  4 

oxp[-g  '21  U  I  as  a  function  of  £,  as  well  as  fe  as  a  function  of  u^. 

It  is  seen  that  the  factor  is  almost  unity  for  l  >  £^.  Figures  7  and  8 
present  the  angular  dependence  of  the  power  spectral  density  for  Z<  0. 
and  l  >  for  various  values  of  the  parameter  a  . 

C .  Covariance  Function  of  the  Ocean  Surface 

The  surface  covariance  function  R(Ar,<))  is  given  by  the  Fourier 
transform  of  the  power  spectrum  of  the  wave  elevation  as  given  by  Kq. 

* 

This  relation  for  is  an  approximate  analytic  expression  for  the  wind 
speed  dependence  of  obtained  from  the  data  in  Pierson's  report. 
Though  could  depend  on  l,  this  dependence  is  ignored. 
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.1 

S  —  cm 


*Th«  run#**  rhown  «nd  th«  v»luoi  of  und  V  »f»  for  u.  *  24  cm/i. 
SOURCE  R«Uf«nc»  1. 


SA-261S-5 


FIGURE  5  POWER  SPECTRAL  DENSITY  OF  OCEAN  SURFACE  WAVE 
ELEVATION  FOR  A  RANGE  OF  WIND  SPEEDS 


(2.58).  Thus, 

m  (  n 

l  SU)< 

'O  ^'-TT 

The  integral  over  $  may  be  performed  analytically  and  expressed  in 


R(£r,g 


,*>  =  / 


o)  exp[.j£&v  cos(«i  -  g)]  d«>  d,£ 


oi  uesaei  tunetions  ot  even  oruer,  .  me  uuiuii.0  of  tne  inicgiui 

’  2m 

are  presented  in  Appendix  A. 


(3. 

terms 
i  ion 


33 


FIGURE  6  PLOTS  OF  exp'l-g2/2C2U4]  AS  A  FUNCTION  OF  C  AND  ^{u.)  .A3  A  FUNCTION  OF 
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In  the  further  development,  we  shall  need  only  the  second  moment,  which 
has  the  form 


"2  “  Ko(u*)  COS"  0  1  °„(U*) 

=  (F  I-  G  )  cos2  9  (  G  sin2  Q  .  (3.17) 

\  o  o/  o 


The  details  of 


the  forms  of  F  (u.) 

o  * 

son  with  Eq.  (2.69) 


the  evaluation  of  the  spectral  moments  as  well  as 

and  G  (u.)  are  also  given  in  Appendix  A.  On  compari- 
o  * 

it  is  evident  that 


and 


0 


(3.18) 


(3.19) 


(3.20) 


Equation  (3.20)  is  the  result  of  the  choice  of  wind  direction  made  in 

the  definition  of  the  power  spectral  density.  The  values  of  H  , 

M.,  as  well  as  G  for  Q  =  0  are  plotted  in  Figure  9  as  a  function  of  fric- 
6  o 

tion  velocity  and  wind  speed.  (Note  that  F  +  G  =  M  for  9  =  <0  Of 

o  o  « 

course,  M  is  independent  of  angle  9,  The  angular  dependence  of  M  is 
o  ^ 

seen  from  Eq.  (3.17),  The  angular  dependence  of  M  and  M  is  complicated; 

6 

however,  they  all  have  maximum  values  at  Q  =  0. 

It  is  of  significance  to  consider  these  numerical  values  for  the 
moments,  along  with  the  Taylor’s  formula  for  R(Ar,0)  given  hy  Eqs,  (2,56) 
and  (2.61). 
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FIGURE  9  WIND  SPEED  DEPENDENCE  OF  SPECTRAL  MOMENTS 


K(Ar,0)  -  Mo  -  i  M2Ar2  I-  £  «,i/  -  ^  «6 


(PoK 


(3.21) 


Consider  the  approximate  representation  oi'  R(Ar,0)  by  only  the  first  two 
terms  of  Eq.  (3021).  Taylor's  remainder  theorem  indicates  that  the 
maximum  error  resulting  from  this  truncation  is  given  by  the  next  term 
after  the  truncation.  For  u,  =  192  cm^s  and  A r  =0.1  cm,  the  terms  in 
the  series  are  numerically  given  by 


12  -3 

-  Mv&r  =  3.64  X  10 


12  -5 

—  M.Ar  =9.63  X  10 
24  4 


16  -5 

120  V*  =3-34x10 

2 

Thus,  M  "  0/2)  M  appears  to  be  a  very  accurate  representation  of 
^  o  fe 

R(Ar,8)  for  displacements  Ar  up  to  0.1  cm. 


D,  Intensity  Distribution  of  Optical  Radiation  Scattered  by  the  Sea 

The  ensemble  average  sputial  coherence  function  is  given  by  Eq. 
(2.30)  as 

(v(Au,Av))  =  exp[-jk(Au  sin  B  cos  A  t  Av  sin  6  sin  A)1 

X  expJ-k2(cos  3  1  cos  B)2[R(0,0)  -  R(Au«4v)jj  .  (3.22) 

It  is  argued  in  Appendix  13  that  the  factor  [R(0;0)  -  R(Au,Av) I  in 

2 

the  exponent  in  Eq.  (3,22)  may  be  replaced  by  (1'2)JI  ir  when  k  is  the 
wave  vector  of  optical  radiation.  Thus 
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<Y(Au,Av)>  =  exp£-j  ~  (^UXD  *'  ^vy0)] 


-k  2  2 

X  exp  — —  (cos  g  i  cos  13)  M(>Ai’ 


(3.23) 


From  Eq.  (3.17), 


("o  "  °o) 


2  2 
Au  l  G  Av 
o 


(3.24) 


From  Eqs.  (2.31),  (3.23),  and  (3.24)  the  average  Intensity  distribu¬ 
tion  of  optical  radiation  scattered  by  the  surface  of  the  ocean  is  given 


2Z  P  cos  a 
o  r 


K)! 


CO 


expj-  ~  (cos  B  F  cos  B)2  ^F  t  G  j 


2  2' 
Au  F  G  Av 
o 


X  exp[-jk(Au£  F  AvT])]  dAudAv 


(3.25) 


where 


x  cos  3  F  :<  cos  B) 
o  / 


(3.26a) 


71  =  (f  +  7J  =  h  (y  cos  0  "  yo  cos  B) 


(3.26b) 


We  note  that  x  and  y  are  the  coordinates  of  the  laser  beam  at  the  tar- 
o  o 

get  plane.  Performing  the  Fourier  transformation,  we  get 
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(3.28) 

This  general  expression  for  the  intensity  distribution,  may  be  specialized 
for  the  backseat  ter  case  by  making  the  substitutions  x  =  xq,  y  =  y  ,  and 
3  =  B. 


(I(x,y))  =  - - 

4rrh“ 


2 

G  x  F 
o 


2  r  2  /  \  2 

Z  P  cos  p  G  x  F  (F  F  G  )y 

o  r _  (  _o _ \  o  o/ 

/Yg  F  F  )g  2h2G  (F  t  G 

\l\  oo/o  L  o\  o  o/ 


(3.29) 


Z  P  cos  p  r  2/2  .21 

.  ,  , .  or  tan  8  /  cos  ff  sm  o 

(Kr/,3)>B  -  2  r  exp  -  2  (  f  f  G  *  G  j 

‘lnh  *\  /fG  f  FlG^  \  o  o  o  ) 


2  2  2 


2,.  2  2, 


(3.30) 


where  the  substitutions  h  /(x  r  y  )  =  tan  3;  x  /(x  t  y  )  =  cos  «,  and 
2  2  2  2 

y  '(x  F  y  )  =  sin  n  as  evident  from  Figure  1  have  been  made.  From 
Figure  1  it  is  also  seen  that  the  same  expression  results  if  the  laser 
transmitter  and  receiver  are  situated  at  the  origin  in  the  receiver 
plane  and  the  position  of  the  illuminated  spot  on  the  ocean  is  described 


by  the  spherical  coordinates  (r^,  a,  n  -  $). 


{•'or  normal  incidence  (3  =  0),  Eq.  (3.30)  leads  to 


(1(0,0))  A  i  (xl  ) 

’  B  -  peak\  */ 

=  Z  P  'Urrh’S/G  (p  *  G  )] 
o  r  [_  y  o\  o  o/ 


(3.31) 


'Vo  note  tliat  I  is  a  function  of  the  wind  speed,  since  F  and  G  are 

peak  o  o 

functions  of  u^.  For  a  moderate  wind  speed  of  7.4  m/s,  corresponding  to 

_ _  _2 

u,  =  24  cm/s,  /G  (F  !  G  )  =  3.47  X  10  .  Assuming  a  range  of  one  km, 

*  o  o  o  ’ 

a  laser  power  of  one  watt,  and  a  reflectivity  (of  the  water)  of  one  per- 

2  -7  2  -2 
cent,  the  factor  (Z  P  /4’-'h  )  has  a  value  of  3  x  10  volts  in  .  Thus, 
or 

the  peak  backseat tered  intensity  under  these  conditions  is  8.66  x  10 
2  -2  -8 

volt  m  ,  or  l,i5  X  10  watts  per  square  meter. 

At  any  other  wind  speed  we  may  write 
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(3.32) 


In  Figure  10,  this  ratio  of  the  peak  intent  j  to  the  peak  intensity 
at  an  arbitrarily  chosen  friction  velocity  (u^  =  24  cm's)  is  plotted 
against  and  U  It  is  noted  that  with  increasing  wind,  or  increasing 

roughness  of  the  surface,  the  peak  intensity  decreases. 

2 

If  we  approximate  cos  3  by  unity  in  Kq.  (3.30),  it  can  be  seen  that 
contours  of  constant  backscuttered  Intensity  are  ellipses  described  by 


2 

cos  g 

V  '  °o) 


2  1 

stn  a  .  2 

~0  3 
O  J 


constant 


(3.33) 


or 


2  2 
x  y 

— *  =  constant 


7F  •  G  \  G 

^  O  Of  o 


(3.341 
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FIGURE  10  WIND  DEPENDENCE  OF  THE  PEAK  INTENSITY  AND  ANGULAR  SPREAD  OF  BACKSCATTERED 
OPTICAL  RADIATION 


'  VBP5r-n!W  <^CTO>A4I  J  J.  ^iA^’AUwA'iBM 


Such  lines  oi'  constant  backseat tered  intensity — or  isophotes — are 

shown  in  Figures  n  to  15  for  various  wind  speeds.  In  each  figure,  the 

isophotes  arc  shown  in  a  polar  plot  against  <.,  for  values  of 

(I/I  )  =  0.8,  0.6,  l/e,  and  0.2.  The  1/e  isophote  is  shown  dotted. 
peaK 

The  value  (I/I  . )  -  1  corresponds  to  the  origin  in  each  case.  The 
peak 

radial  distance  to  the  isophote  at  the  azimuthal  angle  u-  is  proportional 

to  tan  3,  the  tangent  of  the  angle  of  incidence  at  which  the  intensity 

is  the  shown  fraction  of  the  peak  intensity  at  the  same  wind  speed.  The 

absolute  intensities  at  the  isophotes  may  be  compared  by  using  the  ratio 

l I  , (u.)/I  ( 24 )  |  shown  on  each  figure, 

peak  *  peak 

The  semimajor  and  semiminor  axes  of  the  ellipse  of  the  l/e  isophote 

are  characteristic  of  each  wind  speed.  In  Figure  10,  these  values  are 

also  plotted  as  tan  3 1  and  tan  S|  against  wind  speed.  From  these 

‘max  'mm 

figures,  it  is  evident  that  as  the  wind  speed  increases,  the  intensity 
distribution  becomes  broader  as  well  as  lower  in  absolute  value.  This 
is  consistent  with  the  intuitive  notion  that  as  the  wind  increases,  the 
mean  slope  increases,  spreading  the  incident  energy  over  a  larger  angular 

G  C 

spread.  We  note  that  the  tan  3  and  tan  3l  curves  have  approximately 

max  'nun 

the  same  slope.  In  other  words,  the  eccentricity  of  the  ellipses  changes 
only  to  a  small  degree  with  increasing  wind  speed.  This  implies  that 
with  increasing  wind  the  downwind  slope  increases  only  slightly  faster 
than  the  crosswind  slope. 


At  low  wind  speeds,  the  values  of  tan  3  ere  small  even  at  small 

2 

values  of  (I/I  ,).  So  for  these  cases,  the  approximation  cos  3=1 
peak 

is  valid.  However,  at  higher  wind  speeds,  tan  3  becomes  appreciable  even 

for  (I/I  ,)  =  0.8,  and  the  approximation  is  no  longer  valid.  Thus,  in 
peak 

these  cases  the  isophotes  will  bo  ellipses  only  near  the  origin  (small  3)* 
and  at  larger  values  of  tan  3  the  curves  in  Figures  11  to  15  do  not  repre¬ 
sent  isophotes.  The  qualitative  observations  about  the  backseat tered 
intensity  made  above  arc  still  valid. 
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SA-2S18-15 


FIGURE  15  ISOPHOTES  OF  BACKSCATTERED  INTENSITY  AT  u. 


=  1S2  cm/s 


■'.v'^r^-i 


When  cos  0  can  no  longer  be  approximated  by  unity,  and  to  the  ex¬ 
tent  that  the  theory  developed  is  applicable  at  such  angles  of  incidence, 
the  curves  shown  in  Figures  11  to  15  may  still  be  used  to  determine  the 
backscattered  intensity  as  a  function  of  a  and  (3. 

Equation  (3.30)  may  be  rewritten  as 


2  2 
cos  a  I  sin  O' 


2 

tan  3  =  -2  lnj 


I 


2 

I  ,  cos  0, 
peak 


(3.35) 


For  a  constant  value  of  the  right-hand  side,  this  equation  still  describes 
the  ellipses  of  Figures  11  to  15,  with  a  equal  to  the  polar  angle  and 
tan  8  equal  to  the  radial  distance.  In  Figure  16, 


peak 


2 

cos“  3 


]  = 


constant 


Is  plotted  for  the  same  values  of  the  constants  that  are  plotted  in 
Figures  11  to  15.  Tan  0  is  also  shown  as  a  function  of  0.  Thus  at  any 
point  on  the  curves  shown  in  Figures  11  to  15,  we  can  determine  o/  and 
tan  0,  and  from  Figure  16,  the  values  of  0  and  I/Ipeak  mfly  reatl  off. 
It  should  be  emphasised  that  the  validity  of  the  theory  is  questionable 
when  the  small  angle  approximation  is  no  longer  valid. 


E.  Power  Spectral  Density  of  Assumed  Perturbation 

In  order  to  investigate  the  possible  effects  of  various  types  of 
changes  at  the  surface  that  might  occur  as  a  result  of  various  ill- 
defined  and  possibly  unknown  interactions,  the  following  power  spectral 
density  was  assumed  for  a  perturbation  to  be  superimposed  on  the  natural 
"Pierson  spectrum"  discussed  earlier. 


S'a,<5)  &  s'a)F'(o) 


(3.36) 
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SIgiBt5>.-d,--'  -ey^Srp 


F'(0)  =  ~  j\  t-  a'  cos  2(©  -  $)J  for  -n  <  o  <  tt  (3.37) 


I  EGaD  13  /  max 
1  2  eXP)2  (  i' 


s'(j l)  = 
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where 


('  2  ./2\  /  / 2  2 

£  A"  \/8£  f  3£ 
max  y  1 1 _ vo  3 

2  2/1/2  2 

A,  i  n  /\8£  f  3£ 

3  vo  /  \  v  max 


(3.39) 


and  a  =  8.1  x  10  ,  D  =  3.81,  £  =  0.944  cm  ,  and  £  -  8.24  cm"1.  £ 

4  vo  max 

E,  £^,  a^,  and  •>  are  the  variable  parameters  of  the  perturbation  whose 

reference  values  are  chosen  to  be:  £  =  £  ,  E  =  1,  }.'  s  £  .  a'  =  0.5. 

max  3  '  v  vo  1  ’ 

and  0  =  45°.  The  reference  perturbation  is  thus 


S'(£,©)  A  s/(£)E'(0) 
o  o  o 


(3.40) 
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The  following  general  considerations  were  used  to  guide  this  choice. 
The  general  sliape  was  chosen  to  be  essentially  similar  to  the  Pierson 
spectrum.  However,  since  the  capillary  range  (£  =  i  to  10  cm"’1')  is  be¬ 
lieved  to  be  influenced  by  internal  wave  and  other  interactions  with  the 
surface,  the  spectral  peak  of  the  perturbation  was  chosen  to  be  at  the 
beginning  of  the  capillary  region  at  £,^  =  0.944  cm  1,  whereas  the  gravity 
range  is  dominant  in  the  Pierson  spectrum.  In  order  to  obtain  a  few 
finite  initial  moments,  the  spectrum  was  assumed  to  have  a  power  law  be¬ 
havior  similar  to  the  Pierson  spectrum  for  the  viscous  cutoff.  The  spec¬ 
trum  was  chosen  to  coincide  for  E  =  1  at  £  =  l'  =  £  =  8. 2d  cm"1  with 

V  VO 

the  Pierson  spectrum  for  (an  arbitrarily  chosen  but  moderate)  friction 
velocity  of  u  =  24  cm's.  The  angular  dependence  chosen  was  again  sug¬ 
gested  by  Pierson's  report,*  except,  that  the  direction  of  the  maximum 
wave  energy  has  been  chosen  at  an  angle  o  -  •; ;  i.e.,  at  an  angle  *  with 
the  wind  direction.  (See  Figure  3.) 

The  integrated  "energy"  under  the  reference  perturbation  spectrum 

—32  ~32 

is  calculated  to  be  5.804  v  10  cm  ,  composed  of  1.314  X  10  '  cm  in 

-3  2 

the  range  0  to  £  ,  4.  161  /  10  cm  in  the  range  £  to  £  ,  and  2.84 

-5  2  3  v° 

2.84  X  10  cm  in  the  range  £  to  ».  For  comparison,  the  Pierson  spec- 

vo 

trum  for  a  friction  velocity  u^  =  24  cm's  contains  an  "energy"  of 

-32  -52 

8.55  x  10  cm  in  the  range  £  to  £  ,  2.84  x  10  cm  in  the  range 

2  3  V° 

£  ^  to  .o,  and  844.5  cm  over  the  entire  range. 

Let  us  now  consider  the  various  parameters  and  their  effects  on  the 
power  spectrum  of  the  perturbation.  The  defining  equations  have  been 
arranged  so  that  a  change  in  one  parameter  does  not  change  any  ol  the 
others.  The  value  of  E  determines  the  energy  content  of  the  perturba¬ 
tion  with  respect  to  the  reference  perturbation.  The  value  of  £  de¬ 
max 

terminus  the  position  of  the  peak  of  £S'(£)„  }/  determines  the  location 

v 

of  the  cutoff.  ( £  should  always  be  less  than  £7,)  The  range  oi  0  to 
max  v 
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1  in  the  value  of  the  parameter  a'  changes  the  angular  dependence  from 
an  isotropic  behavior  to  a  cosine  behavior  (see  Figure  8).  It  has  no 
effect  on  energy  content.  The  range  of  >j  is  evidently  from  -n/2  to  rr/2. 
However,  due  to  the  symmetry  inherent  in  the  Pierson  spectrum,  it  is  ade¬ 
quate  to  investigate  the  range  0  to  r»  '2. 

We  now  consider  the  perturbed  power  spectral  density 

S  U,o)  4  SO,,«)  i  S'(£,o)  .  (3.43) 

The  perturbed  covariance  function  resulting  from  the  Fourier  transforma¬ 
tion  of  is  evidently 

R  (Ar.9)  =  R(&r,9)  *  R'(Ar,e)  (3.44) 


where 


r"  f  ) 

R^ArjG)  =  J  i  s'  (2)  J  F'(®)  exp[j£Ar  cos(e  -  ®>  I  d©>  d£ 


f  3. 45) 


In  a  manner  analogous  to  the  method  for  Eq.  (3,13)  explained  in  Appendix  A 


CO 

R'(Ar,Q)  =  |  2  s'(£)[joUAr)  -  aj  cos  2(6  -  *)J2(AAr)J  d i 


(3,46) 


Once  again  we  note  that  as  long  as  a'  is  less  than  unity  R7  is  also  mono- 
tonically  decreasing  near  the  origin,  and  the  relative  maxima  do  not 
approach  the  value  at  A r  -  0.  It  is,  however,  interesting  to  note  that 
for  a  value  of  a'  =  1,  the  first  relative  maximum  of  R7  will  occur  at 
Ar  -  6.9  cm.  Whether  R  has  a  maximum  at  this  position  depends  on  the 
relative  values  of  R  and  It'. 
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The  second  moment  of  S'(&,o)  has  the  form  (see  Appendix  A) 

2 

M'  =  F  cos  0  (  211^  sin  9  cos  Q  I  .  (3. "17) 

Using  Eq.  (3.17),  we  may  write 

M,,  *  M'  =  (lfo  *  1^  I  Gq  *  gJ  cos2  e  I-  2H^  sin  9  cos  0 

t  ^G^  (  gJ  sin2  9  .  (3. "18) 

Using  arguments  similar  to  those  employed  in  the  unperturbed  case,  we 
wr  i  te 


<YT<iu,i*>)  -  »*p[-J  ;  (i»xo  ivyj] 

x  exp  j- 1;2 (cos  g  *  cos  B) 2  ^,(0 ,0)  -  I?T(Au,Av)j  j 

■  oxp[-J  n  (4u!to  '  ivyo)j 

Jxp^k2  '2^(cos  g  *  cos  B)2^M9  !  M^)a»'2 J  (3. "19) 


»  exr 


where  we  have  neglected  terms  in  Ar  of  order  higher  than  2  for  reasons 
based  on  the  behavior  of  (y^,)  and  B^,  as  discussed  in  Appendix  B.  It  is 
noted  that  the  expression  for  111^(0,0)  -  RT(Au,Av)j  is  a  quadratic  form 
which  may  be  diagonalized  by  a  suitable  coordinate  transformation.  Using 
Eqs.  (2.31),  (3.49),  and  (3.26),  the  scattered  intensity  distribution  is 
giver,  by 
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<ys,'i>)  =  — 
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xj|exp<~  “(cos  3  i  cos 


s  B)2[(Fo,l,'1lGo,01)Au2i2HiAuAvl  (0oi0>2]| 


Xexp(-jk(§Au*  ?)Av)  1  dAud&v  .  (3.50) 

The  quadratic  form  in  the  exponent  may  be  diagonalized  using  the  rotations 
E7  =  |  cos  ,5  *  t|  sin  {.  ,  An7  =  Au  cos  f  j  Av  sin  5 

(3.51) 

T|/  =  -§  sin  £  '  Ij  cos  5  ,  Av7  -Au  sin  5  »  Av  cos  ? 


to  obtain 


(vs'-’’'’) 


2ZQPr  cos4  0  j  k2 

- - -  exp<-  —  (cos  g  ■  cos  B)‘ 


2f  /  2 

C..u 


x  oxp[-jk(§ 7  Au7  *  n,7  Av7 )  J  dAu7dAv7 


(3.52) 
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(3,56) 


V.'c  note*  that  C,  I),  and  *  are  functions  of  u^,  l  ,  l ' ,  E,  a',  and  ■) . 
Performing  the  Fourier  transformation,  we  get 


Z  P  4 

/  ,_/  „/A  _  o  r _ cos  3  _ 

V'f  *  ’  V  =  2  , _ _  2 

Trh  /cb  (cos  g  *  cos  B) 


X  exp 


L  2(cos  g  *  cos  B) 

For  the  backscatter  case,  it  may  be  shown  that 


2  \  C  D 


(3.57) 


c'  =  2  sin  3  •  cos(o  -  5)  ,  V  =  2  sin  3  •  sin(.v  -  : )  (3.58) 
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(3.59) 


On  comparison  with  Eq.  (3,30)  for  the  intensity  in  the  unperturbed  case, 
it  is  seen  that  the  major  change  is  a  rotation  of  the  distribution  by  an 
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Once  again  we  make  the  approximation  cos  3  =  1  to  rewrite  Eq. 
(3,59) 


2  2 
cos  (o-  -  c )  sin  (if  -  ?) 


D 


tan  a 


cons  tant 


(3.60) 


and,  as  for  the  unperturbed  case,  the  isophotes  are  ellipses  but  tilled 

to  the  coordinate  axes  by  the  angle  *.  (When  the  approximation 
2 

cos  3  -  1  is  no  longer  valid,  considerations  similar  to  those  m  the 
unperturbed  case  apply.)  The  effect  of  the  various  parameters  are  con¬ 
sidered  in  the  next  section. 

It  is  of  interest  to  rewrite  Eq.  (3. '19),  using  Eq.  (3.23),  as 


^.(Au,  Av)  =  <y(Au,Av)>  •  expjj-^k2  2^(cos  g  ‘  cos  B)2  M^Ar“ j 


S  'Y>  W 


(3.61) 


From  the  convolution  theorem,  it  follows  that 

(lT)  =  <1>  *  3{<y'>}  (3.62) 

where  *  denotes  a  convolution  in  the  space  domain. 

From  the  lorm  of  the  functions  y  and  y' ,  both  y  and  y'  are  seen  to 
be  two-dimensional  Gaussians,  even  though  it  is  not  immediately  obvious 
from  the  expression  for  / y' )  since  a  coordinate  rotation  is  required  to 
cast  it  in  the  standard  form.  Their  Fourier  transforms  (lx  and 
then  must  also  be  Gaussian.  The  convolution  of  a  two-dimensional  Gaussian 
with  another  results  in  yet  another  Gaussian  (I  ).  In  this  case,  it  turns 
out  to  be  analytically  more  convenient  to  perform  the  Fourier  transforma¬ 
tion  directly  than  the  convolution. 
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F.  Parametric  Analysis  of  the  Effect  of  Perturbations 

As  indicated  in  the  previous  section,  a  parametrized  power  spectral 
density  was  chosen  for  the  perturbation,  largely  on  heuristic  grounds, 
and  an  analytical  expression  was  developed  for  the  average  backseat tcred 
intensity.  The  variable  parameters  of  the  perturbation  spectrum  are: 

•  i  ,  the  position  of  the  maximum. 

max 

•  i‘ ,  the  position  of  the  short-wavelength  cutoff. 

\ 

•  E,  the  energy  parameter  of  the  perturbation. 

•  a',  the  factor  controlling  the  angular  dependence  of  the 
power  spectrum. 

•  ,  the  angle  the  perturbation  makes  with  the  wind  direction. 

Since  these  parameters  can  be  varied  independently  of  each  other,  it  is 
evident  that  an  infinite  variety  of  conditions  is  possible.  To  limit 
this  variety  to  a  reasonable  number,  it  was  decided  to  choose  a  standard 
condition  and  vary  each  of  the  parameters  about  this  condition.  Table  1 
summarizes  the  conditions  studied. 

In  Figure  17,  is'(jc)  is  plotted  against  log  £  for  various  values  of 

1  and  E  to  show  the  spectral  shape  of  the  perturbations,  for  purposes 
max 

of  comparison,  the  Pierson  spectrum  is  also  shown  for  several  values  01 

u  .  Since  it  is  difficult  to  appreciate  from  this  log-log  plot,  that  the 

*  2  / 
energy  is  held  constant  while  the  parameters  are  changed,  l  S  (£)  is 

plotted  against  log  }  in  Figure  18.  In  this  case  the  area  under  the 

curve,  equal  to  the  energy,  is  the  same  as  the  area  under  the  curve  of 

£S,(£)  versus  £,  even  though  the  spectral  maximum  is  shifted.  Note  that 

a'  and  '•  have  no  influence  on  this  plot  of  the  power  spectrum:  f.e..  for 
1 

all  values  of  a'  and  j ,  the  power  spectrum  of  the  perturbation  is  the 
same  as  the  standard.  Figure  8  conveys  the  effect  f  a'  on  the  spectrum. 
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Tabic  1 

PARAMETERS  OF  THE  PERTURBATION 


Parameter 

Description 

Symbol 

Standard  Value 

Range  of  Variation 

Peak  position 

l 

0.94*1  cm  * 

0.01-10  cm 

Energy  parameter 

max 

E 

1 

0-10 

Short-wave length  cutoff 

i‘ 

V 

8 . 24  cm  * 

1-20  cm  * 

Angular  dependence 

/ 

al 

0.5 

0-2 

Orientation  with 

respect  to  wind 

l) 

45° 

0°-90° 

•)  rotates  this  pattern  about  the  origin  of  the  angle  \).  Note  also  that 
for  E  =  0,  the  perturbation  is  identically  zero,  and  the  results  corre¬ 
spond  to  the  unperturbed  case. 

Figures  19  to  28  show  the  effects  of  changing  the  various  parameters 

on  the  quantities  entailed  in  determining  the  backscattered  intensity 

distribution.  In  Figures  19,  21,  23,  25,  and  27,  I  ,  refers  to  the 

777  7  peak 

backscattered  intensity  at  normal  incidence  (0  =0),  relative  to  the 
value  for  the  standard  perturbation.  Tan  b|  and  tan  gl  respectively 
refer  to  the  semimajor  and  semiminor  axes  of  the  (1/e)  isophote.  In 
other  words,  they  correspond  to  the  value  of  8  when  the  intensity  has 
fallen  to  (1/e)  of  the  value  at  0  =0.  The  angle  §  is  the  angle  by  which 
the  symmetry  axes  of  the  isophotes  are  rotated  from  the  direction  of  the 
wi  nd . 

The  effect  of  the  various  parameters  is  portrayed  in  a  different 

manner  in  Figures  20,  22,  24,  20,  and  28,  where  isophotes  for  which 

.  ....  _  standard . .  .  ... 

i  -  vi'U;  i  ueriveu  l rum  r.q,  u.ou;  are  snown  lor  two  seieeteu 

peak 

values  of  the  parameters,  in  addition  to  the  standard.  The  intensity  on 
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FIGURE  17  POWER  SPECTRAL  DENSITY  OF  THE  PERTURBATION  FOR  VARIOUS  VALUES 
OF  THE  PARAMETERS 


PIERSON  SPECTRUM  FOR  u.  «  24  em/s 


FOR  VARIOUS  VALUES  OF  THE  PARAMETERS 
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FIGURE  23  EFFECT  OF  THE  POSITION  OF  THE  VISCOUS  CUT-OFF 


all  the  isophotos  is  (1/e)  oi‘  the  value  of  the  peak  intensity  (normal 
incidence,  3=0)  when  the  standard  perturbation  is  superimposed  on  the 
Pierson  spectrum  (H  =0).  The  change  in  the  angular  spread  and  the 
rotation  of  the  symmetry  axes  are  to  be  noted. 

While  each  of  these  perturbing  functions  contain  an  almost  insigni¬ 
ficant  amount  of  energy  when  compared  to  the  total  energy  contained  in 
the  Pierson  spectrum  for  u^  =  ?A  cm/s,  it  is  clear  from  these  figures 
that  some  of  these  perturbations  can  greatly  alter  the  scattered  inten¬ 
sity  distribution,  especially  when  they  are  situated  in  the  capillary 
region. 

'or  example,  in  Figure  19  as  the  peak  of  the  perturbing  function  is 
‘a*  id  from  outside  the  capillary  region  to  well  within  the  capillary 
r<-0ion,  the  peak  of  the  scattered  intensity  drops  rapidly,  and  there  is 

i  c‘ 

a  corresponding  broadening  of  the  distribution  as  indicated  by  tan  |j| 

and  tan  b|°  .  Moreover,  while  a  perturbation  outside  the  capillary 

'min 

region  has  little  effect  on  the  symmetry  of  the  scattered  radiation,  a 
perturbation  well  within  the  capillary  region  can  rotate  the  symmetry 
axes  of  the  scattered  radiation  by  a  large  angle  5.  Thus,  it  can  be 
concluded  that,  at  least  for  the  perturbations  assumed  here,  the  scatter!., 
intensity  distribution  is  more  sensitive  to  perturbations  within  the 
capillary  region  than  to  perturbations  at  lower  wave  numbers. 

On  the  other  hand,  Figure  23  indicates  that  moving  the  high  wave- 

number  cutoff  to  even  higher  wave  numbers  has  relatively  little  effect 

on  the  scattered  intensity.  However,  moving  this  cutofl  to  lower  wave 

/  -1 

numbers  docs  affect  the  scattered  intensity,  particularly  for  l  *•-  1  cm 

However,  this  may  be  due  to  the  large  contribution  made  to  the  perturba- 

-1 

lion  by  the  spectral  region  near  1  cm 

Increasing  the  th  (amount  ol  energy)  of  the  standard  perturba¬ 
tion  also  greatly  al  '  • scattered  intensity  distribution,  as 

7‘i 
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indicated  by  Figure  21.  Again  it  is  clear  that  adding  energy  to  the 

capillary  region  oi'  the  spectrum  causes  more  scattering,  as  evidenced  by 

o  o 

the  peak  intensity  falling  and  tan  3  I  and  tan  gl  increasing. 

max  'nun 

We  note  irom  Figure  25  that  changing  the  angular  shape  of  the  per¬ 
turbation  from  isotropic  behavior  (a^  =  0)  to  cosine-square  behavior 
(a'  =  1)  essentially  leaves  the  peak  scattered  intensity  unchanged  but 
affects  the  symmetry  axes  of  the  scattered  intensity  indicated  by  5  and 
makes  the  distribution  more  elliptical. 

Similarly,  Figure  27  shows  that  changing  the  orientation  of  the  per¬ 
turbation  relative  to  the  wind  direction  has  essentially  no  effect  on  the 
peak  scattered  intensity  but  does  affect  the  eccentricity  of  the  ellipti¬ 
cal  in  tensity  distribution  (causing  it  to  be  more  circular)  and  the 
orientation  of  the  symmetry  axes. 

The  change  of  i  with  y  in  Figure  27  is  interesting.  It  shows  a 
maximum  and  decreases  toward  zero.  In  other  words,  even  though  the 
direction  of  the  perturbation  continues  to  greater  angles,  the  intensity 
distribution  returns  to  orientation  along  the  wind.  On  closer  examina¬ 
tion,  it  is  found  that  the  behavior  of  5  at  y  >  '15  is  dependent  on  the 

relative  values  of  F  and  F  in  Eqs.  (.’3. '18)  and  (3. 53)-{3. 56) ,  which  is 

o  l 

governed  by  the  value  of  the  energy  content  of  the  perturbation.  It  is 
found  that  K  returns  to  zero  at  .  -  90  ,  when  E  <.  1.675,  and  that  ;  con¬ 
tinues  on  to  90  at  ;  =  90  ,  when  E  ~>  1.675.  At  E  =  1.675,  V  ~  F^,  It 

i  c  .  e 

should  be  noted,  however,  that  tan  3  and  tan  3  are  approaching  each 

’max  'min 

other,  meaning  that  the  intensity  distribution  is  very  nearly  circular. 

For  this  reason  quite  possibly,  no  great  practical  significance  could  ho 
attached  to  this  observation. 


Comparison  of  these  results  with  those  accompanying  Figures  10-15, 
shows  that,  except  for  symmetry,  the  changes  in  intensity  distribution 


-  -L, ,  ,y  -,yy  . 


brought  about  as  a  result  of  perturbations  are  similar  to  those  brought 
about  by  changes  in  wind  speed  and  therefore  these  two  causes  cannot  be 
distinguished  by  measurements  of  intensity  alone. 

It  has  been  assumed,  of  course,  that  the  ocean  surface  is  homogeneous 
under  the  "scan"  in  botli  cases.  In  other  words,  the  ocean  surface  lias 
been  assumed  to  present  the  same  statistical  behavior  over  the  entire 
scan  area,  which  is  much  larger  than  the  spot  size.  This,  however,  will 
rarely  be  the  case.  Local  fluctuations  will  always  be  present.  Changes 
of  interest  in  the  natural  ocean  surface  due  to  internal  wave,  current, 
or  other  interactions  would  of  necessity  be  inhomogeneous.  While  the 
technique  will  be  insensitive,  unreliable,  or  inapplicable  to  inhomo¬ 
geneities  smaller  than  a  spot  size,  the  variations  on  scanning  over  an 
inhomogeneity  will  be  significant  and  in terpre table  given  appropriate 
additional  environmental  data.  With  the  development  presented  here,  it 
is  difficult  to  assess  these  effects  quantitatively,  lacking  empirical 
data  on  the  spatial  and  temporal  variations  in  the  ideal  spectra  assumed 
in  this  study. 

The  results  >f  Mic-  parametric  study  can  be  summarized  aa  follows. 

The  effect  of  the  perturbation  on  the  intensity  is  most  pronounced  in  the 
capillary  region.  The  energy  content  of  the  perturbation  required  to 
see  significant  changes  Is  on  the  same  order  as  the  energy  content  of  the 
natural  spectrum  tn  the  capillary  region  and  insignificant  compared  to 
the  total  energy  of  the  waves.  The  viscous  cutoff  does  not  have  much 
el  fee1  on  the  intensity.  The  angular  shape  and  the  orientation  of  the 
pertui  .  ‘Ion  al'lect  mainly  the  symmetry  of  the  intensity  pattern  with 
respect  to  the  Aind  direction. 
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IV  SUMMARY  AND  CONCLUSIONS 

The  theory  of  "speckle  patterns”  from  diffusely  scattered  coherent 
light  has  been  examined  and  extended  to  apply  to  scattering  by  the  sur¬ 
face  of  the  sea.  We  have  determined  that  the  surface  characteristics 
can  be  related  to  the  ensemble  average  power  spectral  density  of  the 
scattered  radiation.  However,  since  measurement  of  the  latter  quantity 
is  difficult,  it  would  bo  more  practical  to  measure  the  ensemble  average 
scattered  intensity  distribution.  It  turns  out  that  the  expression  ob¬ 
tained  for  that  distribution  is  similar  to  that  for  high  frequency  radar 
backseat  ter. 

Assuming  a  natural  power  spectral  density  for  the  ocean  surface 
elevation,  the  average  scattered  intensity  distribution  is  determined 
solely  by  the  slope  statistics  of  the  surface  in  the  case  of  optical 
illumination.  The  detailed  distribution  was  calculated  and  its  behavior 
at  different  wind  speeds  was  examined.  We  found  that  as  the  wind  in¬ 
creases,  the  backscattored  intensity  Is  reduced  and  spread  over  larger 
angles.  Wo  also  found  that  the  downwind  slopes  increase  only  slightly 
faster  than  the  crosswind  slopes  with  measuring  wind. 

A  reasonable  parametrized  form  was  assumed  for  the  spectral  shape 
of  a  perturbation,  and  a  parametric  study  was  performed.  The  effect  of 
the  perturbation  on  the  intensity  distribution  is  most  pronounced  in  the 
capillary  region.  The  energy'  content  of  the  perturbation  required  to  see 
significant  changes  is  on  the  same  order  as  the  energy  content  of  the 
natural  spectrum  in  the  capillary  region,  and  is  insignificant  compared 
to  inu  total  energy  oi  the  waves.  the  viscous  cutoif  does  not  have  much 
effect  on  the  intensity.  The  angular  shape  and  the  orientation  of  the 


77 


perturbation  affect  mainly  the  symmetry  of  the  intensity  pattern  with 
respect  to  the  wind  direction. 


Appendix  A 


EVALUATION  OF  SOME  INTEGRALS 


1.  Integration  of  the  Angular  Part  oi  the  Fourier  Transform  of  5(1,0) 


The  Fourier  transform  of  the  power  speetral  density  of  surface  ele¬ 


vation  S(£,©)  gives  the  surface  covariance  function  according  to  Eq, 


(2.58) 


R(Ar)  =|  S(£)  exp(jAr  •  i)  d£ 


(A.  1) 


plane 


Using  Eq,  (3.1),  Eq.  (2.58)  becomes  ir.  polar  coordinates 

f  if  \ 

R(Ar,Q)  =  I  £S(£)  \  I  F(£,o)  exp[jAr£  cos(©  -  0)]  d©>  d£ 


.  (A. 2) 


If  we  use  the  definition 


E  <£)  =  expj^-g2/2i2U<1j  , 


(A. 3) 


Eq.  (3.2)  may  be  written  as 


FU,°)  =  (™^3-*-|l|^l-E(£)f3E(£)aij  jeos  2o+-[l~E(£) Jcos  4cj  .  (A. 4) 


It  is  known1"  that 


<■>■*  [ 


J  (x)  s  -  j  cos(m«  -  x  sin  a)  d*  (m  --  0,1,2,,..)  ,  (A. 5) 

m  r»  i  ■ 


*/*-&* 


r 


V 

b 

5 

i 


From  Eq.  (A. 5)  it  can  bo  shown  that 

tt/2 

2  f 

J,  -  ~  I  cos  2m1?  cos(x  sin  ®)  d<?  (ms  0.1,2....)  .  (A.  6) 

2m  n  I  7  7 

•'o 

By  using  the  symmetry  of  the  trigonometric  functions  involved,  it  can  be 
shown  that 
rt 


I 


m 


cos  2m©  expCjAAr  cos(©  -  9)J  d©  =  (-1)  4  cos  2mQ 


n/2 

■/ 

•'rv 


cos  2m©  cos(£Ar  sin  *)  d© 


=  (-1)  2n  cos  2n;0J  ( ?,&r)  .  (A. 7) 

2m 


Using  Eq.  (A. 6),  and  substituting  in  the  angular  part  of  (A. 2) 


rt 

1 


F(£,©)  expl jAr£  coa(©  -  q)J  d®  =  J^fAr)  -  ^  j'lfl  -  E(i,)]  I  3E(  ^ | 


x  ,I9Uir)  cos  20  -t-  -  [1  -  EU)] 


*  J^(Mr)  cos  40 


(A. 8) 
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*  l^S’i 5^55^S!?9!S5^^^Tji:?^-^i  7&py  jj 


Thus 


R(£r,e) 


CO 

*  J  4S<i)(J„-3 


1  \ 

cos  20.1  I  -  cos  40J  }  d£ 
2  3  4  / 


-05 

i-  |  .eS(je)EU)jj  ^4  «  3a.J  cos  20.J^  -  ~  cos  4© J 


tU 


(A. 9) 


2 .  Evaluation  of  the  Spectr a l  Moments 


th 


The  2n  moment  of  the  power  spectral  density  S (JL,®)  given  by  Eq. 
(3,16)  may  be  rewritten  as 


2n  , 

EU,«)  cos  (0  -  ®)  do>  d£  .  (A. 10) 

IL 

Considering  at  first  the  integration  over  o,  we  note  that 

Tt 


1 


„  _  2n,_  (2n  -  l)::nl2n  „  . 

cos  2m©  cos  (0  -  ©)  d©  =  — -  cos  2mQ  (A.  11) 

2!'(n  -  m) (n  +  m)  ! 

for  0  in  n 


=  0  for  m  >  n 


Substituting  from  Eq.  (A. 4)  and  using  Eq.  (3.18) 
n 


L 


F(A,©)  cos2n{ 


{ On—  I  }  *  *  iTg  4  9  “1 

8-©)d®  =  - - -  -  n(n-l)  cos  y+8r.  cos“  9+2 

2n(n+2):,U  J 

lnE( l)  |a^(n+: 


+2) (2  cos2  e-l)+(n+3) 


8  4  2  1) 

—  (n-1)  cos  0-8  cos  al! 

3  J) 


y  i.  »  i_  eu> 

In  2n 
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(A, 12) 


metna s^u, jtMSi 


2  2  4 

From  Figure  6  we  note  that  E  =  exp(-g  /2i  U  )  may  be  replaced  by  unity 

for  l  >  l  ,  even  lor  the  lowest  wind  speed  (corresponding  to  u,  =12 
1  *inin 

cm/s).  Thus,  for  1  >  .2^  Eq.  (A.  12)  becomes 


n 

1 


2n 

F(.£,o)  cos  (Q  -  o)  do  "• 


i2".--1);',r  (a  cos2  6  .  A  +  „  i] 

2"(„  !■  u:Ll  1  '  J 


I-  -r  IQ  lor  l  >  l 
In  2n  1 


(A. 13) 


Using  Eqs.  (3.3),  which  define  S(£)  over  five  separate  regions,  Eq.  (A. 10) 
may  be  written  with  the  help  of  the  approximation  (A. 13): 

A  t. 


M  =  I 
2n  In 


|  1  vm2n+1  d£  +  i2n  j  1  sL 

***  •'o 

f  s2<mM  d£  +  f 
•4  * 

i 

■ 


U)E(m2nf  l  ax- 


*  (’i„  -  ■  J 


S  (l)/n+1  <14 
o 


,  „  2n+l 

s.UU  d i  i 

4 


l  °° 


. , v  2n+l 

U)l  d t, 


(A. 14) 


3  v 

These  integrations  over  i  may  all  be  performed  in  closed  form  for  n  <  4, 

except  for  the  integral  over  0  -  i  c  which  has  a  series  solution. 

However,  the  series  converges  very  rapidly  and  the  first  few  terms  are 

sufficient.  The  integral  over  £  &  <  »  evidently  does  not  exist  for 

values  of  n  =  4  or  higher.  Thus  only  the  moments  up  to  tr.e  seventh  are 

defined,  and  due  to  the  symmetry  of  H(Ar,f)),  only  the  even  uirmonts  have 

nonzero  values. 
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Adopting  the  notation 


g  =  f  su)i2nW 

"  l  1 


■v  ■  f 1 . 


(£)E(£)^2nhl  d£ 


q  =  s(m2n+1  «u 

n  2 


\ 

i 

'■■I 

Jo 


I  ,  s  2n!-l 
r  =  I  S  (£)*  d£ 
n  I  3 


in  1-1 

s  a  S  ( £)  Jt  d£ 

n  4 


'■  ■  f 

a 

a 

•••I 


t  .  sW^'di 

n  5 


we  have 


(A. 15) 


=  I  g  +  I  h  -i-  /l  +  I  j(q  -f-  r  +  S  +  t  )  .  (A.! 

In  n  2n  n  \  In  2n/\  n  n  n  n/ 


Expressions  for  1,(1,  +  I„  ),  g  ,  h  ,  q  .  r  ,  s  .  and  t  for 

In  In  2n  7  n'  n  n  n‘  n'  n 

n  =  0,  1,  2,  and  3  arc  presented  in  Table  A-l.  The  following  expressions 
for  the  first  two  moments  are  obtained  from  the  table. 


and 


b2  *  [(»  ‘  &'&} 


(A. 20) 


(A. 21) 


On  examination  of  the  terms  in  Eq.  (A. 18)  it  is  evident  that,  if 
a  is  independent  of  },  as  we  have  assumed,  the  terms  are  f unctions  only 
of  wind  speed  or  u^.  Therefore,  wo  rewrite  Eq.  (A. 19)  as 


*2  ’  "„(“*)  5  '  °„ 
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* 


3. 


The  Second  Moment  of  S' (a, ©) 

The  second  moment  of  the  power  spectral  density  of  the  perturbation 


is 


to 

M2  =  j 


' il)  d l  *  j  f'{©)  cos  <9  -  ©)  do  .  (A. 23) 

ro 

Using  Eqs,  <3.37)  and  (A, 11),  we  get 

TT 


1 


2  a  cos  2f  0  a  sin  2i) 

F,(©)  cos  (©  -  Q)  d©  =  - - - -  cos**  Q  +  2  - — 


sin  0  cos  9 


(|-t) 


cos  2<f 


(A. 24) 


Using  Table  A, l,  we  get 


a 

l 


.3  /  .J,  EGAD  11  1 

1  ^  W  di  =  T~  }?  ’  I  CXP  X 


V  +  In  X  + 


I* 

V'  <-x) 

r  •  r; 


\  max'  v  / 


(A, 25) 


where 


X  =  (3/2) 


/i'}2 

\  max  v/ 


(A. 2G> 


Therefore 


» t  t 


I  a^  cos  2i)  ^  I  a'  sin 


M2  «  2 


cos  8+2 


sin  9  cos  9  +  i' 


(f  -  t)  005  2) 


A  * 

S  cos  8  -r  2H^  sin  9  cos  Q  +  g 


(A. 27) 


It  ia  noted  that  F  ,  G  ,  and  II  are  functions  of  l  ,  l* .  Z.  a  .  and  i. 

it  1  max  v  1 
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Appendix  B 


REPRESENTATION  OF  TOE  ENSEMBLE  AVERAGE  SPATIAL 
COHERENCE  FUNCTION 

The  ensemble  average  spatial  coherence  function  is  given  by  Eq. 
(2.30)  as 

(yCAUjAv))  =  expj^-j  |auxq  +  Avy^jJ 

X  expj-k2(cos  8  +  cos  B)2[R(0,0)  -  R(Au,Av)]j  .  (B. 1) 

Consider  the  exponent  of  the  second  exponential  factor;  k  is  the  wave 

2 

vector  of  the  optical  radiation.  For  10.6  nm  radiation  k  has  a  value 

7 

of  3.5  x  10  .  Even  for  angles  B  and  B  as  large  as  60  [when  the  approxi¬ 
mations  made  in  the  derivation  of  Eq.  (2.35)  are  no  longer  valid]  the 
cosine  factor  is  on  the  order  of  unity.  Hence  the  exponential  has  ap¬ 
preciable  value  only  for  very  small  values  of  [R(0,0)  -  R(Au,Av)].  con¬ 
versely,  if  e  is  an  arbitrarily  small  number, 

expj-k2[R(0,Q)  -  R(Au,Av)jj  <  e  (B.2) 


when 


11(0,0)  -  R(Au,Av)  > 


-fin  e/k2) 

\  i 


(B.3) 
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1. 


Behavior  of  (v(Au,Av))  Near  the  Origin 

Using  £q.  (2.61),  the  condition  (B,3)  becomes,  near  fir  =  0, 


12  1  4  /  2\ 

-  M„fir  -  M^fir  +  ...  >  -fin  e/fc  J 


2  2 


24  4 


CB.4) 


If  we  neglect  for  the  moment  terms  in  fir  of  order  higher  than  2, 

^ .  «/2 

fir  >  ( -2  In  s/W„k? 


(-2  In  e/Mgk2^ 


CB.5) 


If  we  choose  the  smallest  value  for  M  (see  Figure  9),  which  will  be  the 


-99 


value  of  G  (0  =  90  )  for  u^.  =  12  cm/s,  and  if  we  choose  e  =  10  , 


fir  >  0.027  cm 
o 


As  discussed  in  Section  II-C,  at  this  small  separation  the  contribu¬ 
tion  of  the  additional  terms  in  fB. 4)  are  negligible,  which  justifies 
disregarding  those  terms  in  this  discussion  and  in  the  expression  for 
[R(0,0)  -  R(fiu,fiv)]. 

Thus  the  spatial  coherence  function  (y(Au>&v))  is  a  very  narrow 
sharply  peaked  function  near  the  origin  and  may  be  represented  accurately 
for  small  fir  by 


<Y(fiu,fiv)>  =  exp j-j  |  ^fiux^  +  Avy^J 


T  2 

Y  CYn  1  ..U  ( 

U 


2  1  21 
nntt  IS  hy*  I 

2  ‘2-  j 


fftsr*  Air*  /'ns- 

(B.6) 
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Behavior  of  (vtAu,Ay))  for  &ry>  -  Ar  <  a 


We  retail  the  Integral  expression  for  JUAr,6)  given  by  Eq,  (A. 9), 
Since  £(£)  £  1  for  all  l,  by  combining  the  two  terms  in  Eq.  (A. 9)  we  may 
write 

CO 

a(Ar,e)  £  J  XS(-j6)[jo  -  aL  cos  2Q  .  J^j  d i  .  (B.7) 

The  value  o i  measured  by  different  investigators  varies  between  0  and 

0.5.  The  variation  of  cos  S  is  between  -fl  and  -1.  Therefore  (J  +  J,) 

o  1 

and  (J  -  J, )  bound  the  variations  of  [J  -  a.  cos  2g  •  JJ,  The  values 
o  1  o  1  2 

of  EJ  (iAr)  ±  J,C£Af)3  are  plotted  as  a  function  of  lAr  in  Figure  B-l, 
o  1 

It  is  seen  that  the  maximum  positive  excursions  occur  for  J  -  J  for 

O  c» 

£Ar  =0,  6.9,  13.3,  and  so  on.  Further,  the  maxima  at  other  than 
Ibx  -  0  do  not  attain  a  value  greater  than  0.6. 

These  observations  mean  the  following  in  interpreting  Eq.  (B.7). 

At  Ar  =  0,  the  Integrand  is  essentially  £S(£),  since  (J  ±  J  )  is  unity 

O  a 

for  all  values  of  t.  Tlius  R<Au,Av>  has  the  maximum  value  at  the  origin. 

As  An  increases,  the  oscillations  in  (J  x  J  )  become  important  and  in- 

O  4 

fluenee  the  value  of  the  integral.  From  Figure  4  it  is  seen  that  S(f) 
is  positive  for  all  values  of  l  and  has  a  fairly  sharp  maximum.  As  Ar 
increases  from  zero,  the  product  iS<A){JQ  i  Jg)  decreases. 

It  is  intuitively  clear  that  the  integral  will  have  a  minimum  (nega¬ 
tive)  value  when  the  first  minimum  of  (J  i  J  )  approximately  coincides 

O 

with  the  maximum  of  the  function  jJS(f).  The  minimum  occurs  at  (iAr)  =  3.5 

for  (J  -  J„)  and  at  (£Ar)  s  5.2  for  (J  +  J„>.  If  the  maximum  of  £S (/,) 
o  2  o  2 

occurs  at  1  ,  tnen  tnc  minx mum  values  oi  the  integral  occur  iior  u  =  l) 

max'  1 

at  Ar  =  (3.5 /l  )  and  (5.2 /f  )— for  the  directions  9=0  and  w/2 

min  1  max  Rax 

respectively.  Similarly,  the  integral  has  a  maximum  value  of 

Ar  =  (6,9 /£  )  and  (8.5 /l  ).  For  instance,  for  u.  -  24  ems, 

max  2  max  max  * 
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FIGURE 


£S<£.)  has  a  ssaxtBiia*  at  *  -  0,0012  css  \  Thus  in  the  direction  6-0, 

the  first  maximum  of  R(fir,g)  occurs  at  approximately  57.5  actors.  Noting 

that  the  oeak  value  of  (J  -  at  6,9  Is  0,6,  it  is  surmised  that  the 

o  2  ' 

value  of  the  integral  will  be  less  than  0-6  R<0,0),  and  condition  (B.3) 
is  amply  satisfied. 

The  conclusion  to  bo  made  from  the  above  discussion  is  that  R(Au,fiv) 
is  a  function  that  smoothly  decreases  away  from  fir  »  0,  whore  it  has  an 
absolute  maximum,  and  whose  further  maxima  are  significantly  less  than 
its  value  at  the  origin.  This  implies  that 

expj-k  lft(0,0)  -  R{fiu,fiv)] j  <  10 

(or  less  than  any  other  arbitrarily  small  value  we  wish  to  impose)  over 
the  entire  (&u,fiv)  plane  cxcopt  near  the  origin,  where  it  may  be  approxi¬ 
mated  by  Eq.  (B.6). 

Since  Eq,  (B.6J  also  becomes  arbitrarily  small  over  the  rest  of  the 
(Au,fiv)  plane,  we  may  extend  the  validity  of  Eq.  (B.8)  over  the  entire 
plane.  Therefore,  when  k  is  large,  we  may  write 

(y(A«,Av)>  a  cxpj»j  ~  (&ux0  *  Avyo)j 

x  oxpj^k^/2  j(cos  g  +  cos  B)2  J  .  (B.B) 

Similar  arguments  may  be  made  in  the  case  of  the  spatial  coherence 
function  of  the  perturbed  ocean  surface. 
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